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^ . Abstract 

' Modified similarity renormalization (MSR) of Hamiltonians is proposed, that per- 

■ formes by means of flow equations the similarity transformation of Hamiltonian in the 
particle number space. This enables to renormalize in the energy space the field theo- 
\ retical Hamiltonian and makes possible to work in a severe trancated Fock space for the 
0^ ■ renormalized Hamiltonian. 

Original works of Wegner in solid state physics has served us as the guiding principle 
in constructing the MSR scheme. 

The renormalized to the second order effective QED^^i Hamiltonian on the light 
. front is obtained. This Hamiltonian reproduces in |ee > sector the standard singlet- 

triplet splitting of positronium and recorves rotational symmetry of the canonical theory. 
The lowest and next-to-lowest energy states of positronium are almost independent on 
the cutoff when in both cases the 'same' (state and sector independent) counterterms are 
5^ i included. 

The electron (photon) mass counterterms are IR (collinear) finite if all diagrams to 
the second order, arising from the flow equations and normal-ordering Hamiltonian, are 
taken into account, and vary with UV cutoff in accordance with 1-loop renormalization 
group equations. 

Both approximations of perturbation theory and Fock space trancation are under the 
control and can be improved systematically within the proposed renormalization scheme. 



*Based on the talk presented at International Workshop "New Nonperturbative Methods and Quantization 
on the Light Cone", Les Houches, France, Feb.24-March 7, 1997 
^E-mail address: elena@hal5000.tphys.uni-heidelberg.de 



1 



1 Introduction 



About two years ago Glazek and Wilson proposed a new renormalization scheme for Hamilto- 
nians, called similarity renormalization. The main idea of their approach is to transform the 
regularized bare Hamiltonian, with a large cutoff A, into an effective Hamiltonian that has no 
matrix elements with energy jumps larger compared to the small cutoff A < A. 

The initial field theoretical Hamiltonian contains usually fluctuations (or states) of all mul- 
tiple energy scales, that couple each other. The effective Hamiltonian of Glazek and Wilson 
has by construction a band-diagonal form in the energy space. This means that the different 
energy scales present in the initial Hamiltonian are decoupled in the effective Hamiltonian. 
This is the main idea of renormalization group approach to take under the control the effects 
of high-energy states, that gives rise to the effective low-energy theory. 

In a relativistic field theory the Hamiltonian usually contains interactions which change the 
particle number of a state. In the Hamiltonian matrix in Fock space representation many-body 
states couple to few-body states. This means that relativistic Hamiltonian is infinite in particle 
number space. The same is true for the renormalized in the energy space effective Hamiltonian 
of Glazek and Wilson as far as it contains the particle number changing interactions. 

What kind of problems arise in solving the relativistic bound state problem by diagonalizing 
finite dimensional Hamiltonian matrices numerically. Suppose, that the energy renormalization 
is performed in the frame of Glazek,Wilson and we are left with the finite in energy space but 
infinite in particle number space effective Hamiltonian. 



First the calculation of effective Hamiltonian so far done are always perturbative [§],|n 
This means that a kind of Tamm-Dancoff (TD) trancation, which restricts the total number of 
particles in any intermediate states, is performed. There are still effective induced interactions, 
together with the canonical one, which change the particle number by two in the band-diagonal 
effective Hamiltonian. To restrict the Fock space in this case is to throw away the interactions 
corresponding to the diagrams with more particles in the intermediate states. Therefore the 
narrow range of the cutoff (band) A, ma^ < A < ma for QED P|, is always introduced in 
Glazek, Wilson aproach, where the correct and cutoff independent values of physical masses 
are obtained. As far as the cutoff A shifted below this range, the perturbative theory blows 
up and the contribution of higher orders in coupling constant diagrams is needed to cure the 
situation [^]. At least for the theories well defined in IR domain as QED we need to elaborate 
the systematic procedure that avoids these difficulties. As was mentiond above the source of 
these difficulties are the interactions which change the number of particles in a state. 

Second, in Glazek,Wilson approach the elimination of far-off diagonal matrix elements is 
performed in each particle number sector, including the sectors that conserve the number of 
particles. In particle number conserving sector the interactions are induced, corresponding to 
the marginal relevant operators of the theory. The elimination of such terms can cause the 
convergency problems 0. Within the framework of Wegner's flow equations, this has already 
been observed for the one-dimensional problem and solved by the proper choice of the generator 
of unitary transformation . This has served us as a guiding principle in constructing the MSR 
scheme. 

The source of all these difficulties is the presence of interactions in the effective Hamiltonian, 
which change the number of particles in a state and mix different Fock components. This 
causes also the general problems of the so-called 'sector-dependent' (occuring in a naive TD 
approximation |12|) and 'state-dependent' (occuring in the method of iterrative resolvents ||15|| ) 
counterterms. Namely, it means,that the counterterms are needed to be added which depend 
on the Fock sector [|12| (and physical state [0) to get the cutoff independent physical masses 
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of states. 

Therefore our aim is to modify the similarity scheme of Glazek, Wilson in a way to de- 
couple different particle number states, that means to eliminate the interactions which change 
the number of particles in a state. Generally speaking, we aim to simulate the effects not 
only of high-energy states on the low-energy one, what conventional renormalization procedure 
does,but also the effects of many-body states on the few-body one by a set of effective interac- 
tions which do not change particle number and do not couple to high-energy states either. In 
this way we are able to control the approximations of the perturbation theory (or any other 
trancation made), and improve it systematically. 

How does the 'modified similarity renormalization' (MSR) work? Say we have the two-body 
bare Hamiltonian, which changes the state of at most two particles (the case of light-front QED, 
QCD). The Hamiltonian is therefore the pentadiagonal block matrix in the representation of 
different number of particle-hole states. Also there is a large number,restricted by the UV-cutoff 
A, of states in the energy representation within each block fig.^ so that there are (infinitely) 
many 'energy states' and 'Fock components' in the field theoretical Hamiltonian. The idea of 
the MSR is to perform similarity transformation, analogous to similarity transformation in the 
'energy space' of Glazek, Wilson, but in the 'particle number space'. This results in effective 
Hamiltonian with a block- diagonal structure in a 'particle number space'. How does it work 
technically? We apply the flow equations in a way to eliminate the far-off diagonal elements in 
the energy space (the matrix elements with large energy jumps \Ei — Ej\ > A) only for those 
blocks that change the number of particles. This is to be done continuously in a differential 
way for all cutoffs from A equal to bare cutoff up to A tending to zero. For the value of A = 
all interactions, canonical and induced, which change the number of particles are completly 
eliminated (except for diagonal in the 'energy space' matrix elements, that do not contribute 
in physical processes). One can say that instead of band-diagonalization in the 'energy 
space' in each 'particle number block' fig.|] (similarity renormalization (SR) of Glazek, Wilson) 
we perform block-diagonalization in 'particle number space', where the elimination of the 
blocks, which change the number of particles, is going in the 'energy space' fig.§ (modified 
similarity renormalization (MSR)). 

Such construction of MSR brings several advantages with respect to the known methods to 
get the effective theories. First, we are able to treat different energy scales in sequence, what is 
usual for the renormalization approach. This solves the problem with divergencies, appearing 
when single unitary transformation (by one step) is performed and all energy scales are treated 



at once [|T2|, [0. Second, in MSR we hope to control both effects of high-energy and many-body 
states. 

One can distiguish then two 'renormalizat ions', in 'energy' and 'particle number' space, that 
are closely related to each other in a complicated way. 

What are the main differences between two approaches (SR and MSR) and consequences for 
the solution of bound state problem? First, in MSR there are still matrix elements with large 
energy jumps in particle number conserving sectors. Second, it is possible in MSR to eliminate 
completly the interactions which change the number of particles, i.e. the limit of A — for 
the effective Hamiltonian is well defined (as A no patalogies occure in MSR in contrast to 
the case of band-diagonal effective Hamiltonian in SR). 

We define the effective Hamiltonian 

i/e//(A) = limi7B(A,A) (1) 

where Hb{X,A) is obtained by the unitary transform of the bare Hamiltonian 

Hb{\, A) = f/(A, A)Hb{A)U+{\, A) (2) 
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The unitary transformation U{X, A) is specified in the flow equation method below. The renor- 
malized effective Hamiltonian is defined 



Hren = hm Heff{A) (3) 

Accepting the above definitions we obtain the renormalized effective Hamiltonian in the 'par- 
ticle number conserving' sectors, which do not depend explicitly on the UV-cutoff. Implicit 
dependence on the cutoff due to the renormalization group running of coupling constants and 
masses is present. We leave aside in this work the question on the new types of countert- 
erms,probably appeared due to the new generated interactions in effective Hamiltonian and 
which can not be absorbed by the running of relevant and marginal operators of the canonical 
theory [0. This means that the physical results become insensitive to the cutoffs. 



The only problem that can arise now by numerical diagonalization of the effective Hamilto- 
nian matrix, that depend on the large UV-cutoff A as size of the matrix and does not depend 



on the artificial cutoff A < A, is the neccesity to introduce the 'sector-dependent' [|T2[ or 'state- 
dependent' 1^ counterterms to get the cutoff independent physical masses of states. This 
problem is solved in the method of MSR, where the resulting effective Hamiltonian is diagonal 
in 'particle number space' and different particle number (Fock) states are completly decoupled. 
Explicitly, the lowest and next-to- lowest energy states are almost independent on the cutoff 
when in both cases the same counterterms are included. This is true despite that we perform 
the calculations perturbatively and work in a severly trancated space [|10| . 



In this work we give the main points of MSR and consider QED on the light front to 
illustrate the method. The key ingredient of MSR are flow equations that perform infinitesimal 
unitary transformation in a continuous way . On the other hand Wegner's flow equations can 
be used in SR scheme to diagonalize a given Hamiltonian in the energy space approximately. If 
in that case one performs the integration of the flow equations only to a finite value of the flow 
parameter, one obtains a Hamiltonian with a band-diagonal structure In MSR scheme we 
use the flow equations to transform the Hamiltonian into a block-diagonal form in the particle 
number space. Both approaches (SR and MSR) are usefuU depending on the problem one wants 
to treat. The purpose of this paper is to illustrate the advantages of MSR manifest in solving 
the bound state problem. 



2 Flow equations and MSR 

The flow equation is written 
dH{l) 



dl 



m.H{i)] (4) 



This is the differential form of the continuous unitary transformation acting on the Hamiltonian 
in general; ri{l) is the generator of the transformation and / is the flow parameter. The aim of 
the Wegner's flow equations is to bring the initial Hamiltonian matrix H{1 = 0), written in the 
energy space, to a diagonal (or block-diagonal where exact diagonalization is not possible) form. 
The finite values of the flow parameter / correspond to intermediate stages of diagonalization 
with the band-diagonal structure of Hamiltonian (SR) ||^. The result of the procedure atl—^oo 
is the 'approximately' diagonal Hamiltonian, where off-diagonal matrix elements do not vanish, 
but they become exponentially small. Also the transformation is designed in a way to avoid 
small energy denominators usually present in the perturbation theory. 
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What is the choice of the generator r]{l) that performs such transformation? Break the 
Hamiltonian into 'diagonal' and 'rest' parts H = Hd + H^. Then the prescription of Wegner is 
i 

^{l) = [HS),HM (5) 

At this step the unitary transformation is defined completly. The only freedom left is the 
princip of separation into 'diagonal' and 'rest' parts. It depends on the problem one wants 
to treat, how to separate the relevant information (degrees of freedom), carried by Hd-, and 
irrelevant one, put in Hr and eliminated by the unitary transformation. 

Our goal is to transform the Hamiltonian into the blocks with the same number of (quasi)particles 
(MSR). This means, that we define the 'diagonal' part Hd as conserving the number of particles 
part before and after interaction, and the 'rest' as particle number changing part. In the 
case of QED(QCD), where the electron-photon (quark-gluon) coupling is present, the number 
of photons (gluons) is conserved in each block of the effective Hamiltonian. 

As a result of unitary transformation the new interactions are induced (see further). They 
correspond to the marginal relevant operators of the theory, i.e. are irrelevant at / = and 
manifest themself only at finite value of /. This means, that these terms must be added to the 
canonical Hamiltonian at / 7^ 0, and therefore give rise to the new terms in the generator of 
transformation ?](/). This in its turn generates new interactions again. 

To be able perform the calculations analytically we proceed further in perturbative frame 
and trancate this series assuming the coupling constant is small. 

As illustration of the method we consider in the work QED on the light front, therefore we 
proceed further in application to it. For any finite value of / one has 

H{1) = Hod + + i/f ^ Jji2) ^ (g) 

where the subscript above denotes the order in bare coupling constant, ~ e"; indices 

'd'-diagonal, 'r'-rest parts correspondingly. The part Hod is the free Hamiltonian, correspond- 
ing to the single particle energies with the structure in secondary quantization a~^a,b'^b{d^d); 
H^^^ denotes electron-photon coupling (of the type a+6+6); Hf^ is the second order diagonal 
part of the Hamiltonian, having the structure b'^d'^bd, b^b^bb, d^d'^dd (in the light front this 
corresponds to the canonical instantaneous (seegull) and new generated to the second order 
interactions in 'diagonal' sector). 

Note, that the diagonal part in MSR scheme is not only the free Hamiltonian but the 
particle number conserving part of the effective Hamiltonian, which contains also corresponding 
interactions. The choice of only Hq for the diagonal part give rise to the band-diagonal structure 
of the effective Hamiltonian in each Fock sector in SR scheme [^. 

The generator of transformation is 

r^(/) = [Hd, Hr] = [Hod, Hi'^] + [Hod, H^^] + ... = 7^(1) + r/(2) + ... (7) 
To the second order the flow equation is written 

^ = [r^, H] = [[Hod, Hi'^]Hod] + [[Hod, H^'W^] + [[Hod, H^'^]Hod] + ... (8) 

The new terms of higher orders e"^ than present in the canonical Hamiltonian are generated by 
flow equations. 

We choose for the basis the single particle states 

Hod\i >= Ei\i > (9) 
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Then in matrix form one has 



Vij = {Ei - Ej)Hll^j + {Ei - Ej)Hlf^j + ... 

^ = -{E. - E^fH^ + [^^'\Hi%, - {E, - E^H^ + ... (10) 

Neglecting the dependence of the single particle energies on the flow parameter Ei{l), one has 
to the first order 

''^"J = -(E-E?H^''^- 
dl 

Hi]^{l) = Hiljil = 0)e-(^'-^^)'' = hI:]^{\ = A ^ oo)e-^^^ (11) 

Here we have used the physical sence of the flow parameter /. Namely, in SR scheme, it defines 
the size of the band A, corresponding to the UV-cutoff, where the effective Hamiltonian is 
defined {\Ei — Ej\ < X) The connection between these two values is 

I = ^ (12) 

In the MSR scheme only the interactions, which change the number of particles, are defined 
in band size A, while the particle number conserving part of the effective Hamiltonian exist 
everywhere (matrix elements with all energy differences are present). 

For the following analyse it is useful to introduce the similarity function 

^•(/) = e-(^-^^)'' = e-^^^ (13) 

which characterizes the behavior (fall off) of the leading order 'rest' interaction with the cutoff 
A. 

We define the initial Hamiltonian at the bare UV-cutoff A ^ cxd (corresponding to the flow 
parameter / = 0) and scale it by flow equations down to the cutoff A < A. The resulting effective 
Hamiltonian is defined in the limit A — (/ — >■ oo), where off-diagonal matrix elements of the 
particle number changing interactions are completly eliminated. It turns out that the diagonal 
elements of these interactions contribute the trivial terms to the 'diagonal' sectors in particle 
space (see further). Therefore one can assume that the effective Hamiltonian at A — has only 
particle number conserving interactions, i.e. has the block-diagonal form in particle space. 

To the second order one has to distinguish the 'diagonal' and 'rest' terms. For the 'rest' 
part one has 

^=[v''\H^%.,-{E.-E,rHS] (14) 

where index 'r' by [r]^^\ H^^^]r defines the particle number changing part of the commutator. 
Introduce 

i^ig(0=e-(^-^^)^'ifig(0 (15) 
i.e. the 'rest' part is defined in the energy band A = 1/Vl. Then the solution reads 

^S(0 = = 0) + /J dU^'~^^^''[v^'\ H^%^,{1') (16) 
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For the 'diagonal' part one has 

^-^=[V''\H^%, (17) 
and the solution is 



^lg(0 = H^d = 0) + jynv^'\H(%,,{i') 



(2) 

Note, that though in general the commutator [[Hod, ]Hod] is not zero, it is not present in 
the flow equation due to the definition of the diagonal part. The corresponding commutator 
[[Hod, H^'^^]Hod\ in the 'nondiagonal' sector insures the band-diagonal form for the 'rest' inter- 
action and also gives rise to the different structure of generated interaction (the integral term) 
in 'rest' and 'diagonal' sectors. 

The commutator [r]^^\ H^^^] gives rise to the new terms to the second order in bare coupling 
e. In the case of QED it induces the new types of interactions, corresponding to the marginal 
relevant operators of the theory, and generates the renormalization group corrections to the 
electron (photon) masses. The coupling constant starts to run to the third order in e. 

Further we consider in 'diagonal' sector all contributions to the second order to protect 
gauge invariance, though in practical calculations we focus separately on the problem of new 
generated interactions and calculation of corresponding counterterms. 

Note, that the MSR scheme (and also SR) enables to build the effective low energy Hamil- 
tonian together with all, 'canonical' and possible to appear 'new', counterterms, found (from 
the coupling coherence condition [T^IjO) order by order in coupling constant e. This defines 



the renormalized Hamiltonian used for the solution of bound state problem numerically. 

3 Renormalized effective Hamiltonian 

Using the flow equations, we derive in this section the renormalized to the second order in bare 
coupling effective Hamiltonian for QED3+1 on the light front. 

3.1 Canonical light-front QED^^i Hamiltonian 

We start with the canonical light-front QED Hamiltonian Hcani devided into free and interacting 
parts 

P- = H,an = J dx-d^x^iHo + Hi) . (19) 

In light-front gauge A'^ = + A'^ = 0, the constrained degrees of freedom A~ and ip- {ijj = 
+ ip-,il'± = A-t?/') can be removed explicitly; this gives the light-front gauge Hamiltonian 
defined through the independent physical fields A_i_ and only ||^ 

Ho = lid^A^WA^ + e (^%T^) ^ ' (20) 



and 

"^667 6^ 



(22) 



7 



H = -ie^ 



'ta-A^^{a-AH) 



(23) 
(24) 



where {o"*} are the standard 2x2 Pauh matrices, and = 29_ = 2^|^. We have used 
the two-component representation for fermion fields introduced by Zhang and Harindranath 
PI '?/'_|_ = {^^. To simphfy the calculations we rewrite all interactions through creation and 
annihilation operators. This turns out to be useful in the fiow equations formalism, f^. 

Following standard quantum field theory procedure we use the momentum-space represen- 
tation for the field operators, and 0, 



= E / ^2(^2^^ (4^..Ae-^- + h.c.) , (25) 

where spinors are Xi/2 — (l^O); X-1/2*'" — (0; 1); with s = —s and polarization vectors e\ = 
^(l,z), = :^(1, — '^); the integration running over the p"^ > only these states, that are 
allowed the light-front theory. 

The corresponding (anti) commutation relations are 

Ka, aJ,A'] = ^9,g'^A,A' , (26) 

where 

V = 2(27r)35(p+ - p'+)5(2)(p^ _ p'^) . (27) 

The light-front vacuum has trivial structure for both boson and fermion sectors, namely aq\0 >= 
0; bp\0 >= 0, simpifying the analitical calculations. The normalization of states is according to 

< Pl,Sl\p2,S2 >= Spj^,P2Ssi,S2 , (28) 

where bp ,.\0 >= \p, s >. 

Making use of the field representation eq. (|25D, we have the following Fourier transformed 

for 

the free Hamiltonian 

the leading order 0(e)-the electron-photon coupling 

Hee^ = E / [9npMil)£xdq+9piP2q{iyxdlq]{bp^hl+K^^^ 

(30) 

where 

^i{Pi,P2,q) = 2— a -a , 31 

q+ pI pI 



The /-dependence of the vertex arises from the dependence of hght-front energies (masses) 
on the flow parameter. Further we have for the instantaneous interactions of the order O(e^) 



n. 



•inst 



S1S2S3S4 "'P1P2P3P4 

^ Xs-j, Xs4, / + _ +\2Xsi Xs2^P3+P4,Pl+P2 

(32) 



and 



inst 

6677 



cr-'cr 



S2 



(p^ - 



vi+<i2,qi+P2 1 



(33) 



here 



~ — g(9+) 
Qjq — ^q.X r~r ) 
V 9"^ 



g(-'?+) 



(34) 



and the 5-simbol stands for the function defined in eq. (|27|), the short notation for the integral 
is 

f _ f dp^d'^p-^ 
Jp^ J 2(27r)3 ■ 



(35) 



In the formulas above we write explicitly the momentum dependence of the coupling con- 
stants as long as / 7^ 0. The initial conditions for the couplings are defined at the value of the 
bare cutoff A — »• cxd(/a = 0), namely 



lim g'^^^ih) = e 

A— >oo 

and for both instantaneous interaction couplings 
hm = ; 

this corresponds to the couplings of the canonical theory. 



(36) 



(37) 



3.2 Flow equations for QED^j^i Hamiltonian on the light front 
3.2.1 Generated interaction in |ee > sector 

Following the procedure outlined in the second section, the leading order generator of the 
unitary transformation is 



r/«(/) 



_ _ Jv^V2(| 



xxt^l{pi,P2, -q)xs^Sg,P2-pi , 



(38) 



VpiPfi^) ~ ^piPfdpiPf 



dgp, 



pf 



^p,Pf 



dl 



(39) 
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where Pi and pf stand for the set of initial and final momenta, respectively, and Ap.^^ = 

J2pT ~ ^Pj ^ ^^'^ light-front fermion energy is p~ — ^ ^ , the photon one q~ — In 
previous notations gp^p^ = Qp^p^q- Further we calculate the bound states of positronium. In 
what follows we consider in |ee > sector 
the generated interaction to the first nonvanishing order 

^eeee ~ 5Z / (0^P3'^P4^P2^piXs3Xs4Xs2Xsi^pi+P2,P3+P4 ) (40) 

with the initial condition limA^oo ^^p"(^a) = 0) 
and the instantaneous interaction 

S1S2S3S4 -^PlP^PSPi 



where 



il-«(^^)=^^ (42) 

The order of the field operators in both interactions satisfies the prescription of standard Feyn- 
mann rules in the |ee > sector. 

To the leading order we neglect the / dependence of light-front energies in the interactions, 
that enables to write the flow equations for the corresponding couplings. 

The flow equations to the first (for the electron-photon coupling) and second (for the in- 
stantaneous and generated interactions) orders are 

^PiPfypiPfK'') 



di 

dVp%{l) 



= (43) 



where 

(44) 

The matrix element < [r]^^\l), Heej] >\ee> is understood as the corresponding commutator 
between the free electron-positron states, namely < P3S3,p4S4\...\piSi,p2S2 >■ 

Note, that the instantaneous and generated interactions are changing with the flow param- 
eter / according to the flow equations in the 'particle number conserving' sector. 

Renormalization group running of both (instantaneous and generated) interactions starts 
to the order 0{e'^), and the electron-photon coupling starts to run to the order O(e^). 

Neglecting the dependence of the light- front energies on the flow parameter I, the solution 
reads 

9piPf{i) = fp^Pf -e + Oie^) 

d(0 = 9;XilA-Q)-e' + 0{e') (45) 



K'.Tfi^) = fdl' < [v^'\n,H,e,{n] >|ee-> +0{e') 
'J 



fpiPf = e =e 
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where the subscript |ee > means, that the commutator is considered in the electron-positron 
sector. The electron-photon interaction exists in the band of size A (|Ap.p^| < A), whereas the 
matrix elements of instantaneous and generated interactions in |ee > sector are defined for all 
energy differences. 

We give further the explicit expressions for the generated interaction, and details of calcu- 
lations can be found in Appendix 0. In what follows we use the notations of this Appendix. 

The matrix elements of the commutator [ri^^\Hee'y] in the exchange and annihilation chan- 
nels are 

^2ii iVpi,P39p4„P2 ~^ VpA,P29pi,P3) 1 

< [7]^^\H^^^] > /Sp^+p^^p.^+p^ = { (46) 

r(an) 



'■^2ii iVpi,-p29p4„-P3 ~'<~ Vp4,,-P39pi,-p2) 



where 

1 dfp^^p^ipj 



Vpi,P2 (0 ^ ■ 



Apjp2 d,l 



^Pi,P2(0 = e-/pi,p2(/) (47) 

and Ap^ p2 = Pi ~ P2 ~ (Pi ^ P2)^ ■ The matrix elements between the corresponding spinors 
in both channels are 

^itf = [xtmPuPs^Pl - P3)Xsi] [Xt^i (-P4, -P2, -{Pl - P3))Xsi] 

(48) 

Mitf = [xt^'l{-Pi,P3, -{Pl + P2))Xs^ [Xt2^i {PU -P2,Pl + P2)Xs^] 

that determine the spin structure of the generated interaction. 

We plug the formulas for commutator [r]^^\ Hee-y] together with the generator r]{l) and cou- 
pling constant g{l), expressed through the similarity function /(/), into the formula for gener- 
ated interaction. This gives rise in both channels 



\Pl P3 ) \ ^Pl,P3 ^P4,P2 



(49) 

poo <^fp^,-p2,^' f J\/ poo 'Vp4,,-V3,^' f J\l 

2i^Aan) -L ( Ja dW JP4,-P3,A'"^ Ja d\' ipi,-P2,A'"^ 



iPl + P2) \ ^Pl-P2 ^P4,-P3 



where in the integral we have neglected the dependence of light-front energies on the cutoff A 
(that is the correction of order 0(e^)), and the connection between flow parameter and cutoff 
(/ = 1/A^) is used. 

Other unitary transformations can be performed, that aim to bring the field theoretical 
Hamiltonian to the block-diagonal form in 'particle number space'. The transformations we 
discuss further, performed in the frame of MSR (section |3.3| ), act also in the energy space and 
differ from the flow equations mainly of how fast the 'particle number changing' interactions 
are eliminated, i.e. by the convergency of the procedure. 

Here we note, that the form of second order generated interaction, induced in MSR and 
written through the similarity functions /p.p^ , is universal for all transformations, which will be 
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discussed. Specifying the similarity function we obtain the exphcit form of generated interac- 
tion, induced by different unitary transformations (section |3.3|) . 
Following the flow equation presciption we use here 



/pi,P2,A = e . (50) 

that gives for the generated interaction in both channels fig.^ 

\/{ex)(\\ _ 2/i^(ex) 1 ^Pi,P3 + ^P4,P2 (T _ f f ) 

^gen *^ ^^Hii / + a 2 1 A2 ' J pi,P3,XJ P4,P2,\) 

yPi - P3 ) ^pi,p3 + ^P4,? 



,P2 



(51) 



^gen K^J — ^^^-^2ii (^+,^+\ /\2 _i_a2 ' V"^ Jpi-p2,XJpA-P3,\l 

\Pl +V2 ) ^Pi,-P2 "I" ^P4 -P3 



By definition, given in the introduction, the renormalized effective Hamiltonian is obtained in 
the limit of cutoff tending to zero (A — 0). In this limit the electron-photon coupling, present in 
generated interaction through the similarity functions fp^pj^, is completly eliminated fpiPfi^ 
0) for Ap.p^ 0, and generated interaction is given by corresponding to 1 in bracket expression 
that does not depend explicitly on the cutoff A. The implicit dependence as renormalization 
group running of coupling constant e and light-front energies (masses) to the next orders is 
present. 

The modified similarity transformation is constructed to avoid divergencies in the form of 
small energy denominators, present in the second order perturbative approach. Namely, one has 
for the generated interaction in the exchange channel in the case of energy conserving process 
i\.Pf = 0) 



^Pl,P3 ~l~ Ap4,P2 _ -Api.PS ~^ ■Ap4,P2 ^ ^ ( I 1 \ /f-r,\ 

A2iA2~Ai2AA 21a A/ 

Pl,P3 ' P4,P2 '-^PiPf ' '^'-^Pl,P3'-^P4,P2 \'-^Pl,P3 ^P4,P2 / 

where 

^PiPf = Pi + P2 ~ P3 ^ P4 ~ ^Pl,P3 ~ Ap4^p2 — Api,-p2 ^ Ap4^_p3 (53) 

due to the total momentum conservation in and 'transversal' directions. The divergent 
contribution in the generated interaction when Ap^ pg ~ Ap^ pj ~ is effectively cancelled by 
the factor in bracket containing similarity functions (1 — /pi,p3,a/p4,p2,a)- 

Also it is obvious from eq. (0), that any energy differences (i.e. VAp-p^) are permitted for 
matrix elements of generated interaction in 'diagonal' sector. 

For the further analyses we write the effective Hamiltonian at the finite value of A, perform- 
ing at the end the limit A — > 0. 

We rewrite the generated to the second order interaction in the form fig.H 



Ai + A2 



A -r 



Af + Ai 



(A) - e N,,, ^4 + ^^4 ■ - 



(54) 
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where we have introduced 



A — Al — Ai . A _ A2 



P1P3 p+ (a;'-x)P+ ' '-'P4P2 p+ (x'-x)P+' 



A _ ^ . A - M'l 



(55) 



(see Appendix ^ for the exphcit definition of these quantities in the hght-front frame). 

The expression eq. (|55|) is written for the rescaled value of the potential, i.e. V\ = Va/P"*"^, 



and the cutoff is defined in units of the total momentum P"*", i.e. A -p+, with / = The 
spin structure of the interaction is carried by the matrix elements M2M, defined in Appendix 
We summarize the instantaneous interaction in both channels to the order O(e^) fig. (^ 



[pT -Pi) 



y{an) _ S - S - 



(56) 



where we have used X^Xi2-'^^siXs4 = SsissSs2S4 + ^siS2^s3S4- For the rescaled potential in the 
light-front frame Appendix 01 eq. (|261D we have 



(57) 



Further we use the generated and instantaneous interactions in |ee > sector, obtained in this 
section, to calculate positronium bound state spectrum. 

3.2.2 Renormalization issues 

As was discussed above the commutator [r]^^^ , -^£67] also contributes to the self-energy term, 
giving rise to the renormalization of fermion and photon masses to the second order. The fiow 
equation for the electron (photon) light-cone energy p~ is 

^ [^7^ ^iHee-yl ^self energy i (^^) 

where the matrix element is calculated between the single electron (photon) states < p' , s'\...\p, s >. 
We drop the finite part and define 6px = p {h)— < \Ho\ >■ Integration over the finite range 
gives 

6p-, - 5p-^ = f' < W'\H^e,] >self energy dl' = JJ^M^^^^ , (59) 



13 



that defines the cutoff dependent self energy 6T,x{p)- The mass correction and wave function 
renormalization constant are given correspondingly, cf. 



d6p 

1 + 



dp~ 



(60) 



The on-mass- shell condition is defined through the mass m in the free Hamiltonian Hq. 

We show further, that to the second order O(e^) the electron and photon masses and 
corresponding wave function renormalization constants in the renormalized Hamiltonian vary 
in accordance with the result of 1-loop renormalization group equations. This can serve as 
evidence for the equivalence of the flow equations and Wilson's renormalization. Therefore 
we have rewritten the mass correction 6m\ through the self energy term, arising in 1-loop 
calculations of ordinary perturbative theory. The negative overall sign stems from our definition 
of the flow parameter, namely for A/ > we are lowering the cutoff dl = — p-dA. 

We start with the bare cutoff mass m\ = -|- 6M^\ where 6Mj^^ is the second order mass 
counterterm. According to eq. (|59D ,eq. (|60|) the electron (photon) mass runs 

ml = ml- [(5S,(m2) - dT.^im')] (61) 

defining, due to renormalizability, the counterterm SMj^^ = 5m\ = — 5SA(m^) and the depen- 
dence of the renormalized mass on the cutoff A 

m\ = m^ + Sm\ = m^ — 5Tj\ . (62) 

We calculate explicitly the self-energy term. The electron energy correction contains several 
terms 



5p-x =<p',s'\H-Ho\p,s>= (^j2^6p-,}j-6'^'\p-p')6ss' • (63) 

The first term is induced by the flow equation in single electron sector, namely comes from the 
commutator [ri{l), Hee'y] 



Spix = - I < [v'^^'K Hee^] >self energy dl' = _^^iAM ; (^64) 
Jlx P 



it reads, cf. eq. ( |280| ) in Appendix^, 

. _ 2 f d^k^dk+ e{k+) ^. + 

xr(p - k,p, -k)r{p,p- k, fe) _ ^_ ^ _ X i'R) ■ (65) 

This term explicitly depends on the cutoff A (/ = 1/A^) through the similarity function, that 
plays the role of a regulator in the loop integration 

i?A = /,%,. = exp|-2(^)'| . (66) 

Eq. (|65|) corresponds to the first diagram in fig. (§). 
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Two instantaneous diagrams, the second and third in fig. contribute the cutofi^ inde- 
pendent (constant) terms. They arise from normal-ordered instantaneous interactions in the 
single electron sector and can be written as 

Spnx = Spn{l = 0) = c„ < 66+ > Vt'\l = 0) (67) 

where n = 2,3 corresponds to the second and third diagrams in fig. (|^), Cn is the symme- 
try factor, < 00+ > stands for the boson {n = 2) and fermion {n = 3) contractions (i.e. 
< a^a^ >= 6{k^)/k^ and < 6p6+ >= 6{p^)), and y^"^*(/ = 0) arises from normal-ordering of 
Hee-yy for = 2 and of -f^eeee for n = 3 (eqs. (|32D and ([33|)). 

These two diagrams 6pn{l = 0) define together with the first one 6pi{l = 0) the initial 
condition for the total energy correction, eq. (p3|). 

Since the diagrams n = 2,3 come from the normal-ordering canonical Hamiltonian at / = 0, 
they must accompany the first diagram for any fiow parameter /. In what follows we use for 
the instantaneous terms the same regulator R, eq. (|66[) 



= ^y^(2;0^^F3IT]X(-^) 

We define the set of coordinates 
X = — 

k = (xp+, xp"*" + /?"*") , (69) 

where p = (p+,p"'") is the external electron momentum. Then the electron self energy diagrams, 
fig. (^), eq. (|284| ) in AppendixjC], contribute 

pUp^, = -^J^x/d^l 



X 



— rri^ ( 2 \ 2m? 

-2 + x\ - — - + 



_4 + /(x) \[x] " ' 7 4 + /(x) ' ' 

f{x) = xm^ — x{l — x)p^ (70) 



and 



oo 



Btt^ Jo 

g2 ^oo 



^TT^ Jo J \ Ix^"^ 



for details we refer to Appendix^. Note, that the transformation in the integrals over x is 
performed before the regulator is taken into account |p. (In the second integral the electron 
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momentum is replaced by the gluon one due to momentum conservation). The brackets '[ ]' 
denote the principle value prescription, defined further in eq. ([77|). 

The loop integral over k eqs. ([70| ) and (|7l|) contains two types of divergencies: UV in the 
transversal coordinate k-^ and IR in the longitudinal component k~^. The physical value of 
mass must be IR-finite. We show, that the three relevant diagrams together in fact give an 
IR-finite value for the renormalized mass; this enables to determine counterterms independent 
of longitudinal momentum. In the wave function renormalization constant, however, the IR- 
singularity is still present. 

Define 

^1 = 1^ ^ (72) 

where P = {P^, P^) is the positronium momentum, p the electron momentum. The transversal 
UV divergency is regularized through the unitary transformation done, i.e. by the regulator R, 
eq. 

Rx = exp |- (^^) I ^ 9{\'5, - |A,,fc|) , (73) 

where the cutoff is rescaled and defined in units of the positronium momentum P+, namely 
A — *• a/^A^/P"*", and Ap^ = p~ — k~ — {p — k)~ = Apfc/p"*". The rude approximation for the 
exponential through a 6'-function changes the numerical coefficient within a few percent; never- 
theless it is useful to estimate the integrals in eqs. ([70|) and ( fflD in this way analitically. From 
eq. (|73|) we have for the sum of intermediate (electron and photon) state momenta (the external 
electron is on- mass-shell p"^ = m?) 

giving for the regulator 

^imax = a;(l - x)X^6i - x^m^ (75) 
and O^nj^^^^) leads to the additional condition for the longitudinal momentum 

^ X ^ Xniax 

Xmax = Z ; 2 I / \ls: \ \ / 

implying that the singularity of the photon longitudinal momentum for x — > 1 is regularized 
by the function Rx. This is the case due to the nonzero fermion mass present in eq. (0) for 
the intermediate state with (1 — x) longitudinal momentum. The IR-singularity when x — > 
is still present; it is treated by the principle value prescription ^ 

11/1 1 , 

(77) 



[k+] 2 \k+ + ieP+ k+ -ieP 
where e = 0+, and P+ is the longitudinal part of the positronium momentum (used here as 



typical momentum in the problem). This defines the bracket '[ ]' in eqs. ([70|) and (^ 
11/1 1 \ 



[x] 2 \x + ij^ X — 



(75 

5i/ 
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Making use of both regularizations for transversal and longitudinal components, we have for 
the first diagram, eq. 




Note, that the third term has the mixing UV and IR divergencies. Combining the three relevant 
diagrams, fig. (^, and integrating with the common regulator, one obtains for the electron 
mass correction 



5ml 
5m\ 



p^{Spi + Sp2 + Sp3 



3m In 



^0) 



The mass correction is IR-finite (that gives rise to IR-finite counterterms) and contains only a 
logarithmic UV-divergency. Namely, when X6i A ^ m 



6m\ 



3e2 2, A2 
— -m In — - 



(81) 

It is remarkable that we reproduce with the cutoff condition of eq. (|7^) the standard result of 
covariant perturbative theory calculations including its global factor 3/8. As was mentioned 
above, the difference in sign, as compared with the 1-loop renormalization group result, comes 
from scaling down from high to low energies in the method of flow equations. 

The similar regularization for the intermediate state momenta in the self-energy integrals, 
called 'global cutoff scheme', was introduced by W. M. Zhang and A. Harindranath [^. In our 
approach the UV-regularization, that defines the concrete form of the regulator R, arises nat- 
urally from the method of flow equations, namely from the unitary transformation performed, 
where the generator of the transformation is chosen as the commutator rj = [Hd,Hr]. Note 
also, that the regulator R, eq. ([T3|), in general is independent of the electron momentum p"*" 
(rescaled cutoff XSi — > A), and therefore is boost invariant. 

For the wave function renormalization constant, eq. (pOf), one has 



d6p 



dp 



577^ Jo 



2 + X x{l — x)2m^ 



that together with the regulator R, eq 



4 + fix) (4 + f{x)y 

results 



x{-R) 



(82) 



In 



X^6, 



2 In 



61 



+ In 



5i 



In 



Si 



+ F In 



X^6i 



X'6i 



In 



A^i (1 XH^ ' 



1 A^^i 



2 X^i + m? 

As X5i ^ A ^ m the function F tends to a constant 



X^di + m? ' A^(5i + m? 
In a; 



dx- 



X — 1 



|A>m 



c 



3 vr^ 
'2^3 



Therefore, by dropping the finite part, we obtain 



1 - 



87r2 



In 



A2 



2 In - + In - 2 - 



M) 



^5) 
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where we have rescaled ^ — > £. The electron wave function renormahzation constant contains 
logarithmic UV and IR divergencies mixed, together with pure logarothmic IR divergencies. 
We mention, that the value of Z2 is not sensitive to the form of regulator applied; the same 
result for Z2 was obtained with another choice of regulator ||^. 

We proceed with renormahzation to the second order in the photon sector. The diagrams 
that contribute to the photon self energy are shown in fig. (^). The commutator [t]*^^-* , ifee7] , 
corresponding to the first diagram, gives rise to (eq. ( P92| ) in Appendix|C|) 



X Tr 



r(fc, k-q, q)Y\k - q, k, -q)] \ r— x {-R) , 



q — k — {q — k:y 



where momenta are given in fig. (0), and the regulator is 



A 



^A = /^,M = exp|-2(^=^j I . (87) 

In full analogy with the electron self energy this also defines the regulator for the second diagram 
with the instantaneous interaction, see fig. (|^), 



2~ 



[q+] J 2(27r)3 ' ' ^ [q+ + k+] 

We define the set of coordinates 

T = X 

q+ 

k = ((I - x)q+ , {1 - x)q^ + K^) 
(q-k) = {xq+,xq^ - K^) , (89) 

where q = {q~^,q^) is the external photon momentum. Then two diagrams contribute (for 
details see AppendixjC], eq. ( |296|) ): 

q^ Sqi = —33^ [ '^^ J '^'^i 
,2 



57?^ Jo 



« (2x=-2x + l) + -j^+(-2 + ^-^)[x(-fl) 



fix) = rri' — x(l — x)(^ 



o2 



1 1 



^ dx \ (i^A ^ (-^) ■ (90) 



^TT^ Jo J \X\ 



Note, that the transformation in the second integral is done before the regularization (by 
regulator the K) is performed 0. 

Making use of the same approximation for the regulator as in the electron sector, we obtain 
for the sum of intermediate (two electron) state momenta 



(91) 



X 1 — X 
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where the photon is put on mass-shell = and the rescaled cutoff A 
used. This condition means for the transversal integration 



Rx 

Xmax 



V2AVP+ has been 



(92) 



and for the longitudinal integration 



Xi < X < X2 

1 — r 

xi -- 



X2 



2 

1 + r 



m 

Wo 



m 



Wo 



(93) 



where the approximate value is used when m <C A. This shows that the condition of eq. (]9ll) 
for two electrons with masses m removes the light-front infrared singularities from x — > and 



X 1. Thus, both UV and IR divergencies are regularized by the regulator R, eq. (p^) 
The mass correction arising from the first diagram, eq. (OT), is 



e2 2 



;a% 



4m' 



2\ 3/2 



Combining together both diagrams with the same regulator, eq. (0), we obtain 



5m\ 



-A 02 r m r — 2m m 



1 



(94) 



(95) 



where r is defined in eq. (pB]). The result shows that the mass correction involves the quadratic 
and logarithmic UV divergencies, i.e. as X62 ^ A ^ m 



5m\ 



-A^ - 2n? In 



The wave function renormalization constant is defined through 



dq~ 



1 /■ , n f 2x2 - 2x + 1 2m'^x{l - x) 
dx I dni < — ^ r— h 



that, with the regulator R, eq. (p2|), results 

/ 2 1 + r 10 8 ^ 

^^ = '-8^(-3^"rr7 + y^+9A^^ 



(4 + f{x)y 



X 



(96) 



-R) , (97) 



(9J 



The photon wave function renormalization constant contains only logarithmic UV divergency, 
indeed as X62 m 



1 - 



3 m?' 



(99) 



and is free of IR divergencies (as is expected from the form of the regulator i?, eq. (PT|)). 
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3.3 Modified similarity renormalization (MSR) by means of differ- 
ent unitary transformations 

We consider in this section the general properties of the unitary transfromation in MSR scheme, 
i.e. unitary transformation performed in the 'energy space' aiming to bring the field theoretical 
Hamiltonian to the block-diagonal form in 'particle number space'. 

Basing on this analyses we construct the general form of the new interactions, generated in 
the effective Hamiltonian to the second order in coupling by such transformation. 

3.3.1 Properties of the unitary transformation in MSR 

In order to illustrate the main idea of MSR consider the flow equations in 'particle number 
space' 1^. Again we break the Hamiltonian into 'diagonal' and 'rest' parts, and choose for the 
'diagonal' part the particle number conserving part of the Hamiltonian H^. and as basis - the 
states with definite number of particles 

H = Hd + Hr 
Hd = He 

Hd\n >= En\n > (100) 

where |n > is the state with particle number 'n'. 
Then the flow equations of Wegner are 

1-1= [Hi, H,] (101) 
One has in the matrix form in particle number space 

dHmn _r rr] _ ( J? ^ J? H 

['/) ^r\mn ^n) ^rmn 

1]inn i^Eyii En) HfYYiYi (102) 

that can be written through the similarity function as follows 

Vmn = jr^{-'^^) (103) 
where the similarity function is defined in 'particle number space' as 

Umn{l) = exp(-A^„/) 

Amn = Em — En (104) 

According to eq. ( |103| ) (due to the second term) the 'rest' part, when m ^ n, behaves with 
the flow parameter / 

Hrmn ^mnHrmn (105) 

and therefore is completly eliminated as / tends to infinity (/ oo, i.e. the band size A — 0). 
One is left with the block-diagonal effective Hamiltonian, where the number of particles is 
conserved in each block, fig. (||). 
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The obvious condition, following from eq. ( |104| ) for the similarity function in 'diagonal' 
sector (m 7^ n) 

Ud = 1 (106) 

is used further. In MSR scheme this means that the similarity function is equal to unity in 
particle number conserving sectors. This condition is equivalent to the choice of the particle 
number conserving part of Hamiltonian for the 'diagonal' part. 

Hd = H,^u, = l (107) 

The idea, that stands behind this excercise, is the similarity transformation of Glazek, Wilson 
but in 'particle number space', which enables to decouple many-body and few-body states in the 
block-diagonal effective Hamiltonian. In the 'diagonal' sectors of the effective Hamiltonian the 
effects of many-body states are simulated on the few-body states by a set of new interactions, 
induced by the similarity transformation and corresponding to the marginal relevant operators 
of the theory. 

The problem, that arise now, is how to simulate also in the effective Hamiltonian the effects 
of high-energy states on the low-energy states. MSR tries to take into account both effects of 
high-energy and many-body states. 

The idea of MSR is again to perform the similarity transformation in the 'particle number 
space' to bring Hamiltonian to a block- diagonal form, with the number of particles conserving 
in each block. Technically, we use the flow equations in the 'energy space', organized in a 
way to eliminate the particle number changing sectors and to generate elective Hamiltonian in 
the particle number conserving sectors. The effects of high-energy and many-body states are 
simulated then by a set of effective interactions, which do not change particle number and do 
not couple to high-energy states either. 

The 'diagonal' sector is the particle number conserving part of the Hamiltonian, and (or) 
the similarity function in the particle number conserving sector is equal to unity 

Hd = H,^ Udij = 1 (108) 

The flow equations are written in the basis of single particle states \i >, i.e. 

Ho\i>=Ei\i> (109) 

where Hq is free (noninteracting) part of Hamiltonian. 

We write the flow equations explicitly in the perturbative frame, i.e. break the Hamiltonian 

i/ = i/od + E(^l"^ + ^i"^) (110) 

n 

where H^"'^ ~ e" and the indices 'd' ('r') correspond to 'diagonal' ('rest') sectors (here we do 
not refer to the definite field theory). The flow equation and generator of transformation to 
the 'n'-th order are 

= Uv^'^Hf^ + H^--'^^] + Y.[[Hit\ Hi--% Hod] + [[Hod, Hod] 
fc k 

r/W = [Hod, Hi--)] + Y.\Hf\ H^""'^] (111) 
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This is written in the matrix form in the basis of states |i > as 



dl Y'' ' 

-{E, - E,) Y.[hP, H(-'%, - (i?, - E,rH^ (112) 

k 

In 'diagonal' and 'rest' sectors one has corresponding 
'diagonal sector' 



f- = Uv^'\Ht'^ + + Y.[[Hf\ Ho,h (113) 



dl 

'rest sector' 

+ E[[^f ' H^''-%^ Hmh - (E. - E^fH!^ (114) 

k 

The main difference between these two sectors is the presence of the third term in the 'rest' 
sector 

-{E:.-E,fHi] = '^.t^ (115) 

that gives rise to the band-diagonal structure in the 'energy space' for the 'rest' part of effective 
Hamiltonian, i.e. 

-^rij ^rij -^rij 

When / — s> oo the 'rest' part is completly eliminated, except for the diagonal in 'energy space' 
matrix elements i = j, which do not contribute to physical values (see further). One assumes 
therefore, that different Fock states are decoupled in the block-diagonal effective Hamiltonian. 

In the 'diagonal' sector the matrix elements with any energy differences are present, Udij = 1. 
According to eq. ( |113|) the new terms are induced in this sector, containing in the lowest Fock 
sectors both the information on high-Fock components (included due to the presence in the 
canonical Hamiltonian the interactions that mix different Fock components) and high-energy 
states. 

The renormalized effective Hamiltonian contains the nonlocal new interactions, correspond- 
ing to marginal relevant operators, and local operators - counterterms, defined order by order 
in coupling constant. One is able therefore to perform at a time renormalization programm 
and work in the lowest Fock sectrors to solve bound state problem. 

The property of the unitary transformation in MSR scheme, we are going to use in the 
next section, is the difference between the 'diagonal' and 'rest' sectors, given by the following 
expressions for similarity functions in both sectors 

Udij = 1 

Urij ^ 1 (117) 
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3.3.2 Generated interaction in MSR 

Here we calculate the interactions in light-front QED^j^i, generated to the second order in 
coupling by different unitary transformations of MSR scheme. Flow equations and similarity 
transformations of Glazek, Wilson, acting in SR scheme, are summarized in Appendix ^ 
The second order interaction, generated by flow equations in MSR in the 'diagonal' sector, 

is 



dl 

-'E,-E, dl ^^^^^ 

where the indices 'd' and 'r' define 'diagonal' (Fock state conserving) and 'rest' (Fock state 
changing) sectors corresponding. The commutator can be written 

[V''\HI%., = - E -f-^ + -f-^^ I (119) 

Neglecting the dependence of energy on the flow parameter, one has the following generated 
interaction 

<(A) = H^iX = A ^ oo) + E + JZ e I (^20) 

where the connection / = 1/A^ is used. Making use of the flrst order solution through similarity 
function 

fij{\ = A ^ oo) 



= = A - oo)^.^^^^^^^^^ (121) 



one has 



Hj^jiX) = H^{\ = A ^ oo) + 5: if;ii(A = A ^ oo)/7« (A = A ^ oo 

k 

h d\' Jjk ja^ ^ Jx d\' Jik \^ J"-^ I ^22) 
Ei — Ek Ej — Ek 



Note, that flrst and second order flow equations eq. (|118|) in 'diagonal' sector, i.e. when the 



condition Udij = 1 is implied, coincide for both MSR eqs. ( |112| ) and ( |113| ) and SR eq. ( P34| )a4 
of Appendix]^ schemes. Assuming Udij = 1 for other similarity transformations of Appendix 

we obtain the same form of interaction, namely eq. (|122[) , generated to the second order in 
'diagonal' sector by any unitary transformation in MSR scheme. 

We proceed along the same line in the light-front frame. Basing on calculations of section 
|3.2.1| , we get the following generated interaction in MSR scheme (for exchange and annihilation 
channels) 



2/i^(ex) -I- I Jx d\' /P4,P2,A'"^ Jx dy /pi,P3,A'"^ 



{Pl -Pa) \ \uP3 ^; 



P4,P2 



(123) 



T/(an)/\\ _ ^2/i^(a») | J\ d\' Jp4,-P3,A'"-^ , Jx dX' Jpi.-p2,A'"^ 

Vgen K.^) - e M^ii 7-+— +T V + V 

\Pl +¥2) \ ^Pl,-P2 ^P4,-P3 
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One obtaines the exphcit form of generated interaction, induced by one of the unitary trans- 
formations in MSR scheme, specifying the similarity function fij in eq. ( |124| ) (Apendix 
Namely, one has for the flow equations of Wegner [|I| 

fpiP2 ~ ^P1P2 ~ ^ ^ (-'-2^) 



and for the similarity transformation of Glazek,Wilson two alternative forms P],[p!^ 

J PlP2 "'piP2^ 

/piP2 ^PlP2 (-'-^^) 

where 

WpiP2 = ^(A-|Ap,pJ) (126) 
We use the general form of generated interaction eq. ( |124D in the light front frame for the 



further analyses in section 4.3 



3.4 Renormalized to the second order effective QED^+i Hamiltonian 
on the Hght front 

In this section we summarize the results obtained in the previous sections. Basing on the flow 
equations we have completed the renormalization of light front QED (LFQED) together with 
calculation of effective Hamiltonian to the second order in coupling constant. 

Consider first the situation when the 'rest' sector is not completly eliminated, i.e. canonical 
electron-photon coupling (the matrix elements) is still present in the energy band \Ei — Ej\ < A. 
This means that perturbative photons are present in the effective Hamitonian and one is able 
therefore to formulate the diagramatic rules for the perturbative expansion in coupling. 

The matrix elements of the renormalized to the second order effective Hamitonian, namely 
the diagrams in different Fock sectors, are depicted in Table |l|, and corresponding analytical 
expressions for 'diagonal' and 'rest' sectors are listed in fig. (^. The diagramatic rules are 
obtained by direct calculation of matrix elements between free particle states. 

Instantaneous diagrams as light front gauge artefact are not included in MSR procedure 
[0],|T1|. This means that instantaneous terms (also belonging to the 'rest' sectors) stay intact 



with the unitary transformation (i.e. do not obtain the similarity factor). 

The matrix elements of the 'rest' (Fock state changing) sectors are sqeezed in the energy 
band A^.p^ = \J2p7 ~ J2Pj\ < -^i the 'diagonal' (Fock state conserving) sectors, exist for 
any energy differences. In the Table ^ we denote as dot the zero to the second order matrix 
elements. 

Using flow equations one has to the second order in coupling the following matrix elements 
in the 'diagonal' sectors between the Fock states 



\ee >—^ \ee >, \eeee >— > \eeee > 

.2 A J xij -^^ dX' ^P4P2A'""> J A dX 



-exIVhij^xO '-rr +T A 1 X 

[Pl - Pi] \ \lP3X \,P2X 

|ee7 >^ |ee7 >, |ee77 >— > |ee77 >, ... 

(roo ^fpikjjj_ r i\i roo ^2^2^ f ^\/^ 
JX dX' Jp2k2X'U-A JpikiX'^^ I (T07\ 
A ^ A ' ^^^'> 
^PikiX ^P2k2X 
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and in the 'rest' sectors between the Fock states 



|ee >^ \eeee >, \eeee >- 

1 



ee >, 



oo 1 



[Pi - Pt. 



^ f 



PiPfX' 



dX' 



roo 1 "•'P4P2 

•1^ U-^.y dX- 



PiPfX' 



A 



V 



P1P3A 



V4P2A 



|ee >^ I77 >, I77 >—y \ee > 

( 1'°° 1 '^fpikiX' f i\l rc 

Jx 7 dX' /p2fc2A'"^ Jx 

exfp.PfxM2ij,xe 



00 1 



df, 



\ 



+ 



fp 



iPf> 



A, 



P2k2>. 



:i28) 



All momenta are given in fig. (|^) We plug the similarity function 



e a2 



^p.Pf = J2Pi -T.P7 



(129) 



into eq. (|127D ,eq. (|128|) to get the explicit form of interactions in both sectors listed in fig. (|]). 
'Rest' diagrams are drawn schematically, to show the difference between the interactions in 
'diagonal' and 'rest' sectors. Namely, for the 'rest' sectors we imply, that the corresponding 
momentum exchange must be done in the diagrams of fig. to get analytical expressions for 
diagrams depicted in Table |l|. 

In the diagramatic rules we write explicitly the dependence of the electron (photon) mass 
on the cutoff 



ml 



rrin 



5S) 



(130) 



where 6J2x is the self energy term to the second order (mo = for a photon); we drop the 
subscript A for the polarization vectors e and spinors x- To the next, third, order one has 
eA = eo(l + 0(e2)). 

The sense of the Table is transparent. Elimination to the second order of canonical 
electron-photon vertex with flow equations gives rise to one-particle and two-particle operators. 
The elimination of the 'rest' sectors to the next orders generate many-particle operators. One 
is able therefore to trancate the effective Hamitonian to the few-particle sectors only in the 
case when the perturbative theory expansion is true. 

As A ^ the 'rest' sectors are completly (here to the order 0(e^)) eliminated and Fock 
states in 'diagonal' sectors are decoupled. One ends up with the renormalized effective LFQED 
Hamiltonian. 

The two-component LF theory, introduced by Zhang and Harindranath [5], as compared 
with four-component formalism of Brodsky and Lepage, is formulated purely in terms of physical 
degrees of freedom; so that each term in the renormalized effective Hamiltonian corresponds 
to a real dynamical process (or give renormalization term). Therefore the 'diagonal' sectors 
of the renormalized effective Hamiltonian contribute for different Fock states to: I7 > - self 
energy photon operator, |ee > - electron-positron bound state (or scattering), I77 > - light-light 
scattering, |ee7 > - Compton scattering, and so on (see Table [|). 

In the next section we use the diagramatic rules given in the Table |1| (at A 7^ 0) to calculate 
the effective electron-positron interaction, which includes generated interaction, instantaneous 
term and perturbative photon exchage. 
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4 Positronium bound state 



4.1 Light front perturbative theory 

The scattering |ee > states are also needed in bound state calculations. Using the propagator 
techniques we include these states where required. We exploit the perturbative theory in the 
coupling constant e, using the diagramatic rules for the renormalized effective theory fig. (^. 

The first order renormalized 667- vertex fpipfH^an contributes to the second order to the 
|ee > interaction term and to the electron (photon) mass renormalization. Physically, it is the 
perturbative photon exchange (photon emission and absorbtion in the case of electron mass 
renormalization), with the energy widths of the photon restricted by the function fpip.x- 



4.1.1 Electron-positron interaction 



According to the light-front Feynman rules the perturbative photon exchange gives rise to the 
following second order |ee > interaction in the exchange channel 



V{1) = gi{l)g2{l)M2^^ 



1 



1 



q^ Pi 



;i3i) 



Pk [-T) Pi -Pk ) 

where Qi stands schematically for the coupling constants in both vertices, namely gp^p2\ = 
efp^p^x and M2ii defines the spin structure of the interaction, coming from the corresponding 
structure of the ee7-vertex; and the two terms in the curly brackets represent two different x"*" 
(time) orderings of the photon exchange with the momentum g, giving rise to the two different 
intermediate states with momenta pk-. Pi corresponds to the initial state. Explicitly one has 
In the exchange channel 



V^t\1) 



-e^{f-P4,-P2{l)xt2^'{-P4,-P2,-q)Xsi){fpi,pi{l)xts^'{PuP3,q)Xs,] 



(132) 



Oipi 



Pt) 



dipt - Pi 



(pi-pt) Pi -Ps -P2 -q (Ps-Pi) Pi -Pi -P4 +q 

with the initial state momentum pi = P = {P^, P^) and momentum transfer q = Pi — Ps, and 

p±2 I jYf2 

^ (133) 



q,pi-p3 



P- 



P+ 



where is the mass of positronium bound state. In the light-front frame holds 

- {P' -Pz -P2 - {Pi-P3Y)d{Pi -pt) = {P' - PY - P4 + (Pi - P^yWiPt - Pi 

A3 



P+(x 



X 



that gives rise for the rescaled potential V = V/P'^'^ 



V, 



PT,X- e^^l^^e^P(^ 



(134) 



(135) 



where A^^i is defined in eq. (|55|) and 



A3 = (kj 



k'±Y + ^{x - x')A + \x 



X 



(136) 
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Here the cutoff A is defined in units of P+, namely A A^/P+. 

Because of tlie absence of Z-graplien in light-front formalism (corresponding to negative 
p^), only one term contribute to the annihilation channel, namely 



(137) 



1 1 " ^ 

X 



where = P and the momentum transfer is q = pi + p2- This gives rise for the rescaled 
potential V = V/P^"^ in the light-front frame, to the expression 

vi-.e=;V.-i,exp{-M±<)j. (,38, 

where and the variables Ai, A2 and Mq , M'q are defined in eq. (^). 
4.1.2 Mass renormalization 

Following light-cone rules the perturbative energy correction of the electron with momentum 
p, coming from the emission and absorption of a photon with momentum k, is 



SPix = J 2{2-k)^ k^^^^^ ~ ^^^3p-k,p,\^\{p - k,p, -k)gp,p-k,xr\{p,p- k, k) 

where 5fee7-coupling constant restricts the energy of the photon. Making use of the explicit form 
for the coupling, one has 

xr\{p - k,p, -k) n(p,p - k, k) ^_ _ ^_ \ X (fl). 

where R = fpp.^ plays the role of regulator. This expression coincide up to the overall sign with 
the energy correction obtained in the previous section from the flow equations method. 

Two instantaneous diagramms, arising from the normal-ordering Hamiltonian, must be 
added to the first term with the same regulator R. Then the full perturbative energy correction 
SPx = SPix + ^P2\ + ^P^x is 

5pl = -5pl (141) 
where is defined in eq. (^). This means for the perturbative mass correction 

5m{^^ = 6Ex (142) 



and the self-energy term is given in eq. (pO[). 

We combine the renormalized to the second order mass, eq. (0), and the perturbative 
correction, eq. ( |142| ), to obtain the total physical mass to the order O(e^) 

ml = ml + 6m^ = (m^ + 5Sa) - = + O(e^) . (143) 
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This means, that to the second order O(e^) the physical electron mass is, up to a finite part, 
equal to the bare electron mass, that stands in the free (canonical) Hamiltonian. 
Along the same line one can do for the photon mass. 

At the end we note, that the similarity function fp^pfX, restricting the electron-photon vertex, 
plays the role of UV (and partialy IR) regulator in the self energy integrals. This means, that 
the regularization prescription of divergent integrals follows from the method of fiow equations 
itself. Moreover, the energy correction (i.e. mass correction and wave function renormalization 
constant), obtained from the fiow equations, coincide up to the overall sign with the 1-loop 
renormalization group result. This is the remarkable result, indicating to the equivalence of 
fiow equations and Wilson's renormalization. 



4.2 Bound state perturbative theory 

In this subsection we define bound state perturbative theory (BSPT). 

First introduce instead of the light front parametrization, used before for the single-particle 
momenta, fig. (^, the instant form 

pif, = {xP+,xP-^ + k-^,Pi) ^ pi^ = {k^,k-^,p'l) = {pi,Ei) 
P2^ = ((1 - x)P+, (1 - x)P^ - k^ ,p^) = (-A;,, -A;^,pO) = (ps, E^) 

Ei = y/p^ + m\ i = l,2 

Ei+kz _ 1 



X 



(144) 

and for the momenta PsjPi the same, but with prime over x,kz,k-^] here x is the light-front 
fraction of the electron momentum, and J{p) is the Jacobian of the transformation: 

^^^^ ^dh^ 2(/+m2)3/2 • 

The instant form system is usefull from the practical point of view: it is easy to see then the 
rotational symmetry, restored in the nonrelativistic case and manifest in positronium spectrum. 

Define now BSPT. We choose the leading order electron-positron potential in such a form 
to simplify positronium bound state calculations. This means, that this potential contributes 
the leading order term to the positronium mass, and perturbative theory with respect to the 
difference between the total second order |ee > interaction, calculated before with the renor- 
malized Hamiltonian to O(e^), and the leading order potential converges. This scheme we call 
BSPT. Surely, our choice is motivated by the form of the renormalized to the second order 
interaction to insure convergence of BSPT. 

In MSR scheme, where the 'diagonal' part is given by the Fock state conserving part of the 
Hamiltonian, the choice of the zero'th order potential for BSPT calculations is quit natural. It 
corresponds to the leading order in coupling of effective Hamiltonian in 'diagonal' sector. This 
insures the convergence of bound state calculations. Note, that in SR scheme the 'diagonal' 
part is given by free Hamiltonian that can not surve for the starting approximation in bound 
state perturbative theory. 

We define the second order renormalized electron-positron potential (e(3)e(4)| Vcoui |e(l)e(2)) 
to the leading order of BSPT in the form of pure perturbative one photon exchange, explicitly 
as the Coulomb interaction 

= - {k^ - k'^y + {k. - k'^y = " (^~p^ • ^''^^ 
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This means that the corresponding leading order Hamilton operator in the |ee > sector is 

H(^) =h + V,oui , (147) 

where h is the free part, defined in eq. (pQl). The wave functions are given as the solution of 
Schrodinger equation 

i/(°)|V^jv(P) >= En\i^N{P) > , (148) 

where P is the positronium momentum, and the eigenvalues and eigenfunctions for the positro- 
nium bound state are defined in standard way of light front frame 

-Cjn — — p+ — 



li^NiP) >= Es,s, ^Jptpt 2(27r)3 5(3)(P - p,p,) $^(x/cxSiS2) (pi) ^^+,(^2)10 > 

^srs2 2(2^)3 ^*Nixk±SiS2) $V(a;/^±siS2) = W (149) 

Mn stands for the leading order mass of positronium. Combining the definitions for the wave 
function and the energy with the Schrodinger equation, we obtain 

^Nix'k'^S.S,) = Y: ^ ,3 V^oul ^Nixk^S.S^) , (150) 



x'{l — X 

or, after change of coordinates according to eq. 

Ml - 4(p' + m") ^m{p's,s,) = Y: I'i V,Up,p') ^n{ps,s,) , (151) 

where the wave function was redefined to have the norm 

Y j d^P^*N{psiS2) ^liPSiS^) = 5nn' ■ (152) 



We aim to obtain the nonrelativistic Schrodinger equation for positronium. Note, that in the 

->2 

nonrelativistic limit ^ « 1 we have 

v^(^W) « i (i - ^±SJ^) 

Mn = {2m + BNf ~ Am'^ + AmBf , (153) 

where we have introduced the leading order binding energy B^^^ . Then to the leading order the 
bound state equation for positronium is 

^ - B„) t.(^.3..) = -Y:Jd'p(^l-^l-V.^,'j 4„(p.,.,). (154) 

Making use of the explicit form for the Coulomb potential, eq. ( |146|) , we obtain the equation 
that determines the leading order bound state wave function: 

P' \ ^ ^ f d?P 



m 



Mp') = ^y (^VpT (155) 
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with 



= $M,s.,«e-(p'S3S4) = <i>f,{p')6s,sjs,s, . (156) 

This is the standard nonrelativistic Schrodinger equation for positronium. Its solution is 
characterized by fi={n,l,m), the usual principal and angular momentum quantum numbers. 
The wave functions are given through the hyperspherical harmonics 

^ Aef ' 

^ ma^ ma ,^ 

B« = = i;r '15^) 

and for the binding energy one has the standard nonrelativistic expression for positronium 
bound state to O(e^). For sake of completeness we write the coordinates used in the solution 

{el = -mBN,p) — > {uo,u) 

Uq = COSCJ = 

u = -p^smu; = „ , ^ , (158) 
\P\ + 

but, for details, refer to 0. 

The electron-positron interaction arising from the renormalized to the O(e^) Hamiltonian is 
given as a sum of two contributions from exchange and annihilation channels V = V^xch + Kmn 
(see explicitly later). We introduce the potential, arising in the nonrelativistic Schrodinger 
equation, eq. ( |155| ), 



JJ{P)J{P') 1 ~ 

V'{p'S'iSi]pSiS2) = hm — {Vexch + Vann)- (159) 

4«i A^T^Y 4m 

Then we define BSPT with respect to the difference 

a I 

dV = V'{p'S3S4;pSlS2) - {-TT^)— ^Ss,ssSs2S4 , (160) 

[P - P ) 



where the leading order contribution is defined in eq. ( |157| ). Note, that, in order to define the 



Coulomb potential, i.e. the ee interaction to the leading order of BSPT, we have taken only 
the first term of nonrelativistic expansion of the Jacobian J{p). 
In what follows we use the matrix elements of 6V, defined as 

< <l>nlJ6V\<!>nim >= J d'pd'p'<^>*^iJp)6V^nlm{P') , (161) 

where ^nim are the Coulomb wave functions given above. 

4.3 Renormalized effective electron-positron interaction 

In this section we consider the properties of effective electron-positron interaction, obtained by 
flow equations and also by the general unitary transfromation in MSR scheme. 
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4.3.1 Renormalized effective electron-positron interaction in light front and in- 
stant frames 

We summarize together all the terms defining the electron-positron interaction, obtained by 
flow equations, in exchange and annihilation channels. 



^exch 



yexch yPT 
yann _^ yPT 



ygen ^ yinst ^ yPT 



ygen ^ yinst ^ y 



FT 
A 5 



(162) 



First we use the light front frame. The generated, instantaneous and perturbative theory 
interactions (rescaled, i.e. V\ = Vx/P^"^) are given corresponding in the exchange channel 



■[/gen 



yinst 



'Ai + 



FT 



(x — x') 
A, 



,\2 ^sisi ^s 



Af+A5 



1 - e ^ 



(163) 



in the annihilation channel 



ygen 



2, A 



inst 



FT 



4e S - S - 

1 <+^'''o 



1 - e 



where in the light-front frame, see eq. ( f^61| ) and following in Appendix |B| 



^ S\S^^ S2Si^\ ' ^2 ^SlS2^SlS3^S2S4,'^^^ 



[X — X 



/\2 



xx'{l — x)(l — x') 



+zmv 2(x' — x) 



Si 



(164) 



SsissSs2S4 ,Ti ■ Sg + 6s-^s.j^6, 
XX 



^2 y± . ^± 



N, 



2, A — '^siS2^S3S4^3 '-^4 ~^ ^ s-iS2^ szSA^ sxss.'^f^ ,/-. n/. , 

xx'[l — x){l — x' 



and 



+im\/2 



S3 



V2 



s,i) 



(165) 



Ti 



\, {k±-k'^y , kl{s2) ^ k'l{s2) 



{x — x') (1 — x) (1 — x') 



kl{s3) , k'liss) 



x' (1 — x') ' 



[l-x) 

'ij 1 



n = 2 



(x — x') X x' 



kiisi) fci(si) 



X 

s = —s 
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with the definitions 

(xk' ~x'ki_? 



Ai 



A, 



+ ni^ix — x')'^ 



A2 = All 



Ai 



Ml 



Ai 



x' — X 



A, 



A, 



X — X 



+ -{x — x')A + |a; — x'\ 



x(l — x) 
A = {kl+m') [j^ 



'2 



X 



f + Ml 
P+ 



1 

X 

P 



k'l + ■m? 
x'{l — x') 

'2 



m 



X 



M, 



N 



2m + B_ 



N 



(166) 



Note, that the rescaled potential, eq. (|162|) , does not depend on the total momentum P"*", i.e. 
is invariant under light-front boosts. 

Second, we use the instant frame. We rewrite the electron-positron interaction in both 
exchange and annihilation channels as follows 



V 



^exch ~l~ '^ai 



+ 



^Ai + A2' 
A? + Ai, 



(1-e- 



+ ^e 
A, 



1 Af+A? 

- \i ^ „PT 



4e2 



,N2 '^■51^*3 ^S2S4 



inst 



[X — X 







+ (4e' 5s, 
where one has in the instant frame 



:i-e- 



+ 



'0 

IS3 '^52 54 





inst 



Ml 



(167) 



X 



1 + 



k. 



X 



1 + 



M' 



(168) 



a/p^ + m? J 
,Q _ 4(/ m^) ; Mq^ = 4(]y^ m' 

and for all the quantities, defined in eqs. ( |165| ) and ( |166| ), the substitution x{kz), x'{k'^) is to be 
done. The symbols cf^" and others are introduced to indicate the origin of the different terms 
(here generated interaction coming from the exchange channel); all c = 1. 

The expression eq. (|167|) has the physical interpretation. First notice, that Ai and A2 de- 
scibe the energy differences (energy transfer) in two corresponding ee7-vertices appearing in the 
electron-positron interaction. Then the generated interaction (cf^*^) contributes mainly hard 
photon exchanges ^ ~ ^ >> 1, while the term arising from perturbative theory c^J' gives 
rise to soft photons. Though the effective electron-positron interaction generally describes 
low energy physics, namely the renormalized effective Hamiltonian is constructed to eliminate 
perturbative (relativistic) photon exchanges, the information on the high energy sector is accu- 
mulated in the generated interaction. This makes possible to interpolate between two sectors, 
i.e. the sum of both terms in eq. ( |167D recovers the whole range of photon energies. 

In the next section we analyse the general properties of the effective electron-positron in- 
teraction, generated by the unitary transformation in MSR scheme. 
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4.3.2 Electron-positron interaction in MSR 

We give here the general analyses of the electron-positron interaction, obtained in MSR scheme. 
Basing on this analyses, we calculate in the next section positronium fine structure. 

We summarize the interactions in |ee > sector, generated by the unitary transformation of 
MSR (section p.3.2|) , in exchange channel fig. (H) 

= -^"^^^^1, + 

^3 

Vr' = -(^^^--^--' (169) 

where the similarity functions 

/ia = /a(Ai),/2a = /a(A2), (170) 

are specified in eqs. ( [I24| ) and ( |125| ) section ^.3.2| ; the energy denominators Ai,i = 1,2,3 
together with A^i ^ are given in eqs. (|165|) and ( |166|) . Introduce the following notations 



Ix dX' 

Then the generated interaction is written 

where we have used the identities 

01 + 02 = 1-/ia/2A, (173) 

The explicit form of similarity functions together with 6i factors is 
for Wegners fiow equations 

/a(Ai) = e-^ 

^1 = Ti^(l-/A(Ai)/A(A2)), (174) 

for the similarity transformation of Glazek,Wilson 

/a(AO = ^^(A^-IAil) 

9i = 9{\Ai\-\A2\)e{\Ai\-\') = 9{Ai-A2)9{\Ai\-\^), (175) 



where 



Ai = ^,A2 = ^, (176) 
x— X x' — X 
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We stress once more, that the unitary transformation of MSR is constructed in a way to avoid 
small energy denominators. Namely, the situation Ai = A2 = 0, corresponding to the process 
with no momentum transfer (Ai = A2 = means in light front coordinates x = x', k± = k'j_, 
i.e. g = pi — = 0), is excluded from the generated interaction for any finite A > 0. Further 
we assume for all positive A 

\//^"(g = 0)=0,VA>0, (177) 

This means also, that the diagonal matrix elements, left after the limit A ^ is performed 
in the 'rest' sector to the first order in coupling {eej vertex when Ai = 0, A2 = 0), do not 
contribute to the second order generated interaction in 'diagonal' sector. Therefore one can 
assume, that the 'rest' sector is completly eliminated by the unitary transformation of MSR. 
In this limit (A 0) the perturbative part of electron-positron interaction is zero 

V^^ = 0, (178) 



and electron-positron interaction is defined by the instantaneous interaction eq. (|169|) and the 
generated term at A ^ 

= -eiV,.|^ + ^ 

VX^O = - ^jf^,/ s.s:A.s, , (179) 

The expression 

y\ee> ^ ygen^ ^ yu^st ^ (^gQ) 

gives the effective electron-positron interaction obtained in MSR. Note, that when A ^ no 
difficulties (compared with SR) on the convergency of unitary transformation (MSR scheme) 
and (or) convergency of perturbative theory occure. This is closely related with the maintenance 
of the rotational invariance and gauge invariance of perturbative calculations. 

The explicit form of generated interaction when A — is given in the case of 

flow equations 

similarity transformation 

V.r„ ^ -e^-N. (^i^ + . (182) 

We consider now several special cases to analyse the general expression eq. ( |169| ). 
I. Consider the energy conserving process (EC), i.e. when the following condition is 
valid 

for the light front coordinates of the process fig. (^). This gives rise to 
Ai = A2 

/a(Ai) = /,(A2), (184) 
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Making use of the formula eq. ( |172| ), one has for the generated interaction 



Vr\EC = + i)(l - /ia/2a) = -e^iVi,,i(l - , (185) 

and for the perturbative part 

Vr\EC = -e'N,,,^fl , (186) 

where 

fix = /a(Ai) (187) 
When the following condition holds 

M^ = ^(Mo' + M;2), (188) 
one has for the energy conserving process (Ai = A2) 
Ai + A2 



2 

This gives rise to 



A3 , (189) 



VriEC = -eXAJ-(l-/A'(A3)) 

^3 

Vx^^Iec = -e'N,,,^fl{As), (190) 

^3 

where we have introduced 

A3 = -P- , (191) 

X — X 

Expression eq. (|190|) shows clearly the sense of the unitary transformation in MSR scheme: 
how the elimination of ee'-f vertex in the 'rest' sector with the flow parameter / (or cutoff A) 
generates the new interaction for the |ee > Fock component in 'diagonal' sector. For the energy 
conserving process together with the condition eq. (|188|) the resulting generated interaction (at 
A — > 0) is given by the pure one-photon exchange diagram (i.e. by the initial perturbative 
interaction at A = A — > cx)) 

VrolEC = VCI^Jec . (192) 

II. We consider the collinear limit, i.e. when |x— a;'| — * 0, and show, that though the instan- 
taneous term is singular in this limit the whole electron-positron interaction is finite. Namely 
the generated to the second order interaction insures the absence of collinear divergencies. 

For this purpose we rewrite the electron-positron interaction eq. ( |169|) as follows 

2 A. 1 

9 /I. 1 1 X 1 \ . 

-^^^H2^a; + a;)-a;]^^-^-^-) 

-eXA^(^-i)(^i-^^2), (193) 



2^Ai A 



2 
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The first two terms (perturbative photon exchange and instantaneous interaction) define the 
initial value of electron-positron interaction (A = A ^ oo), the next terms are generated by 
the unitary transformation and describe the change of electron-positron interaction with A. 

In the collinear limit |x — x'| the A-dependent terms (namely similarity functions 
fix = /a(^), = and 6i factors, i = 1,2) are given at the value A = 0. 

Since the singular contribution is carried by the instantaneous interaction we consider the 
whole interaction only in diagonal spin sector, i.e. 6sj^s^6s2S4,- 

The first two terms can be written 

1 46^ 

2e^ 1 ~ 2 
= -^{7^Ni,x{siS2 S3S4) + A3- w^Ss.s^Ss^s^) , (194) 

A3 ^ {XX) 

where A^'i^a, defined in eq. ( |165[ ), is given 

l.r / N ^{k±-k'^y 1 /,2/ 1 ,'2/ 1 1- 



+ f{SiS2 S3S4) 



/(++ - ++) = ,n VI 

xx'[l — x){l — x') 



f{+-^+-) = kl— - + 



'2 



xil — x) x'il — x) 



\xx' (1 — x){l — x'] 



Q-K't'-'t'') 



k± = fc_L (cos 0, sin {/)) , (195) 

We have used the notations '-|-'-spin up, '—'-spin down. 
Assuming the condition 

M^ = 1(Mo2 + M;2), (196) 

one has 

2 Jk^-k'^? , 1 ,.2, 1 \ ,'2, 1 



As, -T = 2^ ^ + 



(x — x')^ {x ~ x'y {x — x') 



tM(-, 7--)) 

1 — X X 1 — X X 



~^ ^ xx' ~^ {1 — x){l — x')^ ' ^ ^ 



Combining all terms together, we get 

1 4e^_ 



-e^A^i,A(siS2 5354)-^ - - — ^-^—^Ss^si^s 

A3 [x — x' 



2eV 2/ 1 1 N n\ 

A3 \ xx' [l — x){l—x') J 



(198) 



This part is finite, i.e. collinear singularity in the instantaneous interaction is cancelled exactly 
by the perturbative part. To the leading order in 6x = x — x' one has 

1 4e^ 

-e^Ni x{siS2 8384)^ Tv?^siS3Ss2S4 = 

A3 [x — x')^ 

^""^ ^m2(l + —^—,) + fl=AsiS2 ^ 5354)! + 0{5x) , (199) 



{k^-k'^Y\ (l-x)2 
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For the next term we perform the expansion with respect to 6x 

1 {k^ + k'^f 6x^ 



1, 1 1 
+ — 



then one has 

1 

2e2 



+ 0{5x' 



(200) 



1 J_ 



(1 - /ia/2a) 



{k± - k'^Y 



x^fl — x) 



0{5x) 



(201) 



Both expressions eq. (|199[ ),eq. ( p01| ) do not depend on the exphcit form of unitary transforma- 
tion performed. The last term does depend in colhnear hmit on the exphcit form of similarity 
function. 

In the collinear limit ^j-factors eq. ( p.74| ),eq. (|1 75|) are given in the case of 
flow equations 



^il 



A? 



Af + Ai 



and similarity transformation 



^il 



^(Ai - A2) 



(202) 



(203) 



Performing the expansion in 5x, one has in both cases corresponding 



2 VAi A2 



h-e2 



2e^ 



flow equations 



2 , 



(204) 



0, similarity transformation 



Combining all terms, one has the following electron-positron interaction in collinear limit (con- 
sidered only the spin sector Ss^sa^s^gJ 



Vx'^^lx = x') = < 



mH\ 



[i-xY 



+ 



flow equations 



2e^ 



1 



x^(l—x)'' 



+ f\x=x'{SiS2 S3S4) 



+ f\x=x'{.SiS2 S3S4) 



(205) 



similarity trans f or mation 



that is finite. 

III. Let us analyse what kind of interaction (repulsive or attractive) arise from the effec- 
tive Hamiltonian in the electron-positron sector eq. ( p.69| ). 

First note, that in the light-front frame the generated interaction has definite sign for any 
value of A. Namely the term in bracket of generated interaction eq. ( |169| ) is always positive, 
since energy denominators Ai, A2 are positive on the light front, and > 0, /jj > {fij{X) 
has the same behavior (A), defined in Appendix ^). 

When A = generated interaction together with instantaneous term give rise to the attrac- 
tive electron-positron interaction in the whole parameter space. 
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To the leading order of nonrelativistic approximation 

\v\ 

— < 1 , (206) 

m 

one has in the instant coordinate frame 
Ai ~ A2 ~ A3 = A = - p'Y 

[p - p r 
[p - p r 

A = , (207) 

This gives for the electron-positron interaction in the whole nonrelativistic range of A 
X <ti m 

V\ee> ^ -e^ \ - ^^^5 5 , (208) 



the A- independent result. Making use of the following expressions 

^^1 ~ ^ / i\o "sisals 



(fc± - k'^Y 



- p'f = (k^ - k^y + (k, - k,f ^ {k^ - k^f + Am\x - x'f , (209) 
one get to the leading order of nonrelativistic approximation the 3-d Coulomb interaction 



e2 



^l-"> ^ 16m' [-J^Zpy) '^-^35...4 ; (210) 

hence the rotational invariance is restored to this order. This result eq. ( |210| ) is valid for any 
nonrelativistic value of A, and for any unitary transformation in MSR (i.e. V/jj) performed. 

One can interprete the independence of the interaction eq. ( p08| ) on A, when A is in nonrel- 
ativistic domain A -C m. In this region all photons, that give rise to perturbative (relativistic) 
corrections, are almost eliminated. Therefore effectively only the lowest Fock state (|ee >) con- 
tribute to the interaction and trancation of the whole Fock space to the lowest Fock component 
is valid in this case. 

The explicit A-dependence (through the similarity functions) in the interaction signals to the 
presence of the terms in effective Hamiltonian that mix different Fock component, i.e. particle 
number changing interactions. Then trancation is not possible. 

One way out is to use different approximations, resulting in freezing effectively the contri- 
bution of particle number changing (Fock state changing) interactions (nonrelativistic approx- 
imation 1^; the definite choice of scale (size of the band in SR scheme) A, depending on the 
system considered |TT|). 

In the case of MSR the effective Hamiltonian is defined in the limit A ^ 0, where all Fock 
states are completly decoupled and one is able to solve the bound state problem separately in 
each 'diagonal' sector. 
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In the next section, making use of the nonrelativistic approximation for the electron (positron) 
momentum 

b1 



m 



0(a) < 1 , (211) 



in the effective electron-positron interaction eq. (|167| ), obtained by flow equations, we calculate 
the positronium splitting analytically. 



4.4 Positronium's fine structure and rotational invariance 



In the nonrelativistic approximation we obtain the following effective electron-positron inter- 
action, eq. ( |167| ), 



1 1 1 



,1-^1 

2(27r)3 4m2mV 2m'^ J 



1 _ e"^^^^ 1 cff - 



-2( -^Y „PT 



{p — p'Y 

r r inst 

2; 



{p-p'f + \x- x'\{M^ - Ml) 



+ 



1 -e 



-2 



prp 



Ml 



, n ■ - ^ ^ 

where the energy denominators and exponential factors were simplifled using 



J,nst 
-"an 1 



X — X 



2m 



1 + 



-f2 

2^ 



m 



Ai = As = [p-p'f + 0\ m 
A3 



2 ( P_ 

m 



{p- p'f + \x- x\{Ml -Ml) + 0[m^ (- 

m 



Ai = A2 



2m{p' — pY 
{k'z - K) 



M'^ = Am^ + 0\ m^[ — 



1 + 

2\ 



A 



2m{p' — pY 
(k'z - kz) 



Ml = 4m2 + AtuBn + O ( — 



m 



4m2 + 4m5S^ , 



and the explicit expression of Jacobian for the coordinate change is 



1 - 



p' 



2m2 



O 



(212) 



(213) 



(214) 



having introduced the leading order binding energy bP. Making use of the nonrelativistic 
approximation B'^/m « 1 we have for the interaction 
111. ^2 



2(27r)3 4m 2m 



\ 2my 



39 



X 



^gen PT 



1 + 



inst r r 



+ 



S3S4 



(0) 



{P-P') 



|A 



2 PT ^ 



m 



(215) 



where (c^*^", c^'^) showes that both term, generated and perturbative interactions, contribute 
to corresponding term (we remember that all c = 1). 

The remarkable feature of the part of interaction standing in the first bracket is that it 
does not depend on the UV cutoff A. The next term in the second bracket arises from the 
perturbative photon exchange and has the typical 'energy shell' structure for the relativistic 
effects, namely these terms are important when A >> m. Further we calculate the ground state 
positronium mass and therefore restrict the cutoff to be in the nonrelativistic domain 



A « m 



(216) 



where the second term in eq. ( |215| ) vanishes and we are left with the following form for the 
renormalized ee interaction in the nonrelativistic approximation: 



V'{p,p') 



1 1 



2(27r)3 4m 2m 



P \ 



{p-p') 



^ 2m2 j 



(217) 



16e^m^ 



-r> 



'2S4 



+ 



^^^2 , eff 

4^2 ^^a.ri 5 "-an J 



+ 4e2cr*5,,,.5, 



We perform the nonrelativistic expansion of the factors Ni and N2 appearing in the interaction. 
The term A^^i contributes in V' to the order 



2 / -2 

16m2% ( 1 + 1^ ) + ( k\K + A;',* A;; 



-m{si + S2) [Ji^kA] - A{k^ + }^^_f + AsiS2q\ 



imV2{x x) ^^^^ ej^ ■ dg^ss ^S2S4 + _ ^^^^ _ ^/^ ' ^siS2 ^s^s^^ 
= 8 6s^s3 ^S2S4 m {iql - siql) ^1 - 
- 554^2 6s^s3 m {iql - S2q±) 1 + 



m J 
kz + k' 



m 



qz {ipl - SipI) + -S2qz{q± - isiql) 
qz {iPl - S2PI) - ^siqziql - ^s^ql) 
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,/\2 



2m' 



{x — x') 
xx'{l — a;)(l — x') 



While the term N2 contributes to V to the order 

, 2^ 



o fey 



2m' 



xx'{l — x)(l — x') 



32m' 



In these formulas we have used [k'j_, /cj 
following variables have been introduced 



£ij3 and el = --^{s,i)] 



Q± 
Pi. 



kj_ — kj_ 
k± + k\ 



X, y) , qz = k'z- kz 



(218) 
also the 

(219) 



We leave aside for the future work the analysis of the expressions for A^i and N2, where also in 
this form some terms can be identified as spin-orbit and spin-spin interactions in the transverse 
plane and in longitudinal (z) direction. 

Instead we follow Q, where an analogous calculation of singlet-triplet ground state mass 
splitting of positronium was performed in the similarity scheme. This means, that we can, 
except for the leading order term 0(1), drop in A^^i the part diagonal in spin space. Also the 
terms of the type / 
splitting, since 



k^^kz , k^^kz , k^k\^ do not contribute to the ground state mass 



/ 



(220) 



averaging over directions, gives zero. 

We obtain for the ee-potential to the leading order 0(1) of nonrelativistic expansion 



1 1 



2(27r)3 4m2m 
16e^m^ g| 



g2 



1z 



1 + 



p2 

2m? ^ 

P \ 



m^ 



(gen PT\ 



16e^m^ 

1z 



1 + 



pz 
m^ 



jinst 



a 



27r2g2 
V{t) ( 



1 + 



1 



f ^ 

2m^ J 
2m2 ) 



(221) 



Remembering q = p' — P-, Fourier transformation to the coordinate space with respect to q has 
been performed in the last expression. To the leading order of NR expansion we have reproduced 
the Coulomb potential, defined before as the leading order of BSPT. Note, this is true for any 
UV cutoff within the nonrelativistic range A <3C m. 

We combine this expression with the kinetic term from the Schrodinger equation, eq. (|155|) , 
and write it in the form 



1 

m 



1 + 



V{r) 
2m 



+ V{r) . 



(222) 
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Here the potential V{r) plays a different role in the two terms. In the first term, corresponding 
to kinetic energy, it generates an effective mass of the electron, which depends on the relative 
position and manifests the non-locality of the interaction. The second term is the usual potential 
energy, in our case, the Coulomb interaction. 

The energy of the Coulomb level with quantum numbers {nlm) is standard 

=< <l>nlm\V'^'^\<!>nlm >= J pd' p' <^l,^{p) V'^'^ <!>r.lm{P') = , (223) 

where the Coulomb wave functions ^nim were defined in eq. (|157|) . We have used in eq. ( |223|) 
the following representation 

[p - p') = — [u-u) 



{u — u'Y ^ n 



£p = dn, (224) 

and also orthogonality of the hyperspherical harmonics 

j dnY;Y^. = 5,,, . (225) 

More details can be found in p. 
The next to leading order O 

1 1 1 / e2\ 



2{2-KfAm2m\ J 

X {8m{iql - Siql)Ss^s3Ss2S4 - 8m{iql - S2ql)Ss^s3Ss2S4) (226) 



contributes (because of the spin structure) to the second order of BSPT: 

' ~ h W^i ■ ^ ^ 

The order O ^(^) ^ (cf. remark after eq. (^T^)) is 

1 1 1 / \ 
'^^^^^ " 2(27r)3 4m 2m [^^^ + Se^5,^,^5s.,s^5s^s3 + ^e^^s.-s^ds^sA (228) 

and contributes to the first order of BSPT: 

5M2=<$ioo|5V^(')|$ioo> . (229) 
Both contributions were calculated in with the result 

SM"^ = 5Ml + 5Ml 
5 

< ll^Af^ll > = ma^ 

' ' 12 

<2|5M2|2> = <3|5Af2|3 >=< 4|5M2|4 >= -maS (230) 

6 
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where the eigenvectors in spin space are defined as follows: 



> , |4 >= 1+ +> . (231) 

Making use of the relation between Coulomb energy units and Ryd = |ma^ we have the 
standard result for the singlet-triplet mass splitting for positronium, |a^Ryd. The degeneracy 
of the triplet ground state n = 1 refiects the rotational invariance, manifest in the system in 
nonrelativistic approximation. 

5 Conclusions and outlook 

We proposed in this work the new Hamiltonian approach, the modified similarity renormal- 
ization (MSR) of Hamiltonians. By means of fiow equations, acting in the energy space, the 
unitary transformation is performed, which aims to bring the field theoretical Hamiltonian to 
the block-diagonal form in Fock space. 

The renormalized effective Hamiltonian, constructed in such a way, solves two problems. 
First, all 'diagonal' Fock sectors, i.e. the sectors belonging to Fock state conserving, are com- 
pletly decoupled in the effective Hamiltonian. Second, elimination with the flow equations 
of 'rest' sectors, i.e. Fock state changing, perturbative order by order in coupling constant, 
enables to flnd to these orders all counterterms, corresponding to canonical operators (rele- 
vant and marginal) of the initial theory and to new operators, generated by flow equations. 
The 'new' counterterms are combinations of the 'canonical' one or can carry also new types of 
divergencies. This question we stayed aside in this work. 

Summarizing, one is able to answer two questions: the physics of what Fock state and of 
what energy scale is descibed by the renormalized effective Hamiltonian. 

Namely, one can pick out any block from 'diagonal' sector of effective Hamiltonian, acting 
on the deflnite Fock state, and solve separately for this block (Fock state) the corresponding 
physical problem. In particular, it is possible to trancate the block-diagonal effective Hamilto- 
nian to the lowest Fock state (here discussed |ee >), plug the corresponding interaction (from 
'ee' block) into the Schrodinger equation, put it on the computer and solve the bound state 
problem (for positronium bound state) numericaly. 

The problem that arises by numerical diagonalization of the lowest Fock sector matrix in 
the energy space (or by numerical solution of integral bound state equation, see main text) 
is the dependence of physical mass spectrum on the size of the matrix (or on UV cutoff im- 
posed on the energies in the integral). This dependence is to be absorbed by including the 
corresponding counterterms, calculated analyticaly on the previous step. MSR scheme insures, 
that counterterms to add do not depend on the state (i.e. are the same for the lowest and 
next-to-lowest exited states) and on the Fock sector (i.e. trancation to the |ee > state do not 
prevent counterterms corresponding to the diagrams with more particles in intermediate state) 
[|T0|] . This is due to the new interactions and corresponding counterterms, generated by the 
flow equations in MSR scheme. 

In the work we have outlined the strategy to build the renormalized effective Hamiltonian, 
that can be fulfllled by means of flow equations in a systematic way (to the end A ^ 0) without 
knowing apriore any properties of the system considered further (i.e. no approximations are 
needed to make in the procedure). This is the main advantage of MSR scheme. 



II >-^(l+ ->-!-+» 

|2>=-^(|+ ->+|- +>) 



13 >= 
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We remind, that through MSR scheme one simultaneously renormalizes in the energy space 
the initial field theoretical Hamitonian and constructs the effective Hamitonian, for which the 
Fock space trancation is valid. But rather the analytical solution of flow equations, as of any 
other iterrative method, demands to apply perturbative expansion in coupling constant. This 
means, that applicability of perturbative theory is closely related to the possibility of working 
in a trancated Fock space. Within the proposed approach one is able to improve sytematically 
this approximation. 
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A Similarity renormalization 



In this appendix we skip the idea of similarity renormahzation (SR) proposed by Glazek and 
Wilson 0, the renormalization scheme of Hamiltonians by means of unitary transformation, 
and consider two alternative methods: the flow equations (FE) of Wegner [|^ and similarity 
unitary transformation (SUT) that are used in SR scheme and give rise to the renormalized 
Hamiltonian. 

The similarity transformation aims to bring the field theoretical Hamiltonian to the most 
diagonal form, namely only the matrix elements between the free states with 

\Aij\ < A (232) 

where Ajj = Ei — Ej, are present in the renormalized Hamiltonian. For this purposes the 
continuous unitary transformation is performed to preserve unchanged the spectrum (eigen- 
values) of the initial bare cutoff Hamiltonian. The demand of diagonal structure does not 
define completly the generator of the transformation. This freedom is used to eliminate small 
energy denominators in the final renormalized Hamiltonian. This results in a system of two 
self-consistent non-linear differential equations for the Hamiltonian H{1) and the generator of 
the transformation ri{l). The dependence on the continuous flow parameter / in the flow equa- 
tions by Wegner is replaced by the cutoff dependence A in the similarity unitary transformation, 
with the connection 

/ = l/A^ (233) 

in the renormalized Hamiltonians. 

The difference of the two methods (FE) and (SUT) consists of the residual freedom in the 
choice of the direction of the infinitesimal rotation, actually defining how fast the non-diagonal 
matrix elements vanish. 

We summarize the equations for both methods, written in matrix form. 

I. The flow equations by Wegner [||]: 
with 

Uij = exp(-/Aj) (236) 

and 

fij = Uij . (237) 
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II. The similarity unitary transformation by Glazek and Wilson 0: 

^-".^'•l^'h^r.,^^, (238) 

and 

fij = Uij exp(rjj) . (240) 
Also the following transformation is used 0: 

§^^..,1,,//,]., + ^^. (241) 

and 

h = u,, , (243) 
where for both similarity transformations 

u,, = e{X-\A,,\) (244) 

and 

Uij + Tij = 1 . (245) 

Also other choices for the similarity function Uij with the step behaviour are possible 0]. 

Remember, the function f^j defines the solution for the leading order interaction term. The 
first order equations for H and rj, written through the /-function are unique for both methods 
(I. and II.) 

i^ffl(0 = ^S(/ = 0)-^j|§^, (246) 



with the connection given in eq. (|233|) in the renormalized values, and dl — > dX implied. The 
eqs. ( |246| ) and ( p47| ) will be exploited further for the calculations in the main text. 
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B Calculation of the commutator [r]^^\l), H^^^] in the electron- 
positron sector 

Here we calculate the commutator [ri^^\l)^Hee'y\ in the electron-positron sector. The leading 
order generator rf^^^ is: 



xx^3rKPi,P3,-g)Xsi5g,_(p,_p3) , (248) 



where 



Vpipa (0 ~ ^PiP3 ■ 9piP3 



1 



1P3 



^piP3 dl 



(249) 



A 



PlPa — Pi ~ Ps ~ (Pi ~ Ps) ! the electron-photon coupling 

H..^ = i9*p,pS^yx'dq' + fep4(04'«ig') ibpJ>P2 + bpJ-P2 + d-pA2 + d-pJtp,) 

(250) 

where 



rz Pi,P2,g =2— + a - a 



Pi 



(251) 



and the tilde-fields are defined in eq. (0). Further we use the identities for the polarisation 
vectors and spinors 



XsXs' 



A 



Using the commutation relations, eq. (pB]), and identities eq. ( p52| ) we have 

■ 2 I 'ipiP3yp2P4, p+ -p+ 'ipip-iyp2Pi ^+ _ ^+ I 



e67J 



(252) 



(253) 



^. y U,p,^U_p^U-p^Up^ Up^Ll_p^Ll^pgUp2 T fpgU_p^U_p_jfp2 T U^^LL _p^U-p,j^Up^ ) . 

x(X^,rKPl>P3,Pl -P3)X.i) (xiri(p2,P4,P2 -P4)XS2) ^Pi+P2,P3+P4 > 

where the first two terms of the field operators contribute to the exchange channel, and the 
next two to the annihilation channel. We take into account both s- and t-channel terms to 
calculate the bound states. The : : stand for the normal ordering of the fermion operators and 
(|) is the symmetry factor. The sum over helicities Sj and the 3-dimensional integration over 
momenta pj, i = 1, ..4, according to eq. (p5|) is implied. We rewrite for both channels 



iVp 



e(-(p+-P+)) 
^ -ipt-pt) 



P4,-P2 



V-pi-pidpi, 



pay 



iV-p4.,-P29pi,p3 ^Pl,p-i9-pA, 



-P2J 



^S^iS h+ rl+ rl - h 

/\U "pi+p2,p3+p4 "p3Sa P454 P2S2"P1S1 



667J 



(254) 



-M. 



2ij V 2 y 



(an)^l^^ J e(pj+pp 



(Pi+Pj) 



{Vpi,-P29-PA,P3 V-pA,p39pi-p2, 



e(-(p++P+)) . 

+ -(P++P+) l^- 



Pi,Pz9p^~P2 



n 



Pl,~P29-P4,P3y 



h+ rl+ rl - h 

Pl+P2,P3+P4 "p3S3"'p454 P2«2"PlSl 
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where 



(ex) 



2ij 
(an) 



2ij 



{xtXl{-P4,P3, -{Pl + P2))XsJ ixt.^iiPl, ~P2,Pl + P2)Xsi 



(255) 



The first term in the exchange channel with pi > corresponds to the hght-front time 
ordering xf < with the intermediate state = P3 + {pi — Ps)" + P2 ^ the second term 
Pi < Pt s^nd xf > xf has the intermediate state = Pi — {pi — Ps)" + Pa - Both terms can 
be viewed as the retarded photon exchange. The same does hold for the annihilation channel. 

Consider only real couplings and take into account the symmetry 



V- 



■PA-P2 



'4,P2 5 



-PA-P2 



9p4,P2 ■ 



(256) 



Then ^3S3,p4S4| [r]^^\ Hee-y] \piSi,p2S2) , the matrix element of the commutator between the free 
states of positronium in the exchange and annihilation channel, reads 



+P2 ,P3+P4 



-^2n i^Pl^P39p4,P2 "I" Vp4,P29pi,P3) 

n/f an 1 



iVpi,-P29p4,~P3 ~^ Vp4-P39pi,-P2/ 

We rewrite this expression through the corresponding /-functions 



(257) 



Vpi,P39p4,P2 "I" Vp4,P29pi,P3 



l^Pl,P3 



Vp 



1,-P2 



9p4,-P3 "I" Vp4,-P39pi, 



-P2 



dl 



/p4,P2(0 ~^ 



^fp4,P2 (0 



A, 



P4,P2 



A 



P1-P2 



dl 



/p4,-P3(0 "I" 



dl 

1 dfp^^.p^il) 



A, 



P4,-P3 



dl 



(258 

fpi,-P2iO 



with Apj^p2 = Pi ^ P2 ~ {Pi ~ P2)~ ■ As we have mentioned in Appendix A this form in terms of 
the /-function is universal for all unitary transformations. We exploit further this expression 
by specifying the /-function to compare the effective interactions in different renormalization 
schemes (see Appendix C). 

We calculate the matrix elements M2ii, eq. (188), for both channels. Here we follow the 
notations introduced in 

We make use of the identities 



6'^ + ise 



5'^ + ise'^ 



with s = —s and xtXs' = ^ss'] also of 



xt(^'x-s 



-V2, 



(259) 



(260) 



with e* = —Es and = —6^-,. 

We use the standard light-front frame, fig. (|^), 

= (2;P+, xP^ + k±) , P2 = ((1 - x)P+, (1 - x)P^ - 

P3 = {x'P+, x'P^ + k'^) , Pi = ((1 - X )P+, (1 - x)P^ - k'^) 

where P = {P^, P^) is the positronium momentum. 



(261) 
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Then, to calculate the matrix element M2ii in the exchange channel, we find 



P^[xt3^'{Pl,P3,Pl-P3)Xsi] = X. 



S3 



(k± — k\ Y a ■ k\ a ■ k\ x — x 

: V W T^cr + a' + im - 

[x — x') x' X xx' 



Xsi 

(262) 



and 



P^[xtXi-P^^ -P2, -(P4 - P2))Xs,] 



S2 



\^ {k^-k'J ^ (a-k, 



X — X' 



l-x 



(J + (J 



a ■ k\\ X — X 

— im- — rC" 

1 — x' 



{l-x){l-x') 



where we have introduced 
Ti = 



Tl5s,s, + ^"^ (1 -X){1- X') ^~^)^^^^2^^2f4 



Xs4 

(263) 



Tl = 2 



X — x' {1 — x) (1 — x') 
(A;^ - k'J' kiisi) A;'i(si) 



(264) 



X — 



X 



X' 



and 

Finaly we result 



(265) 



{x — x') 



,/\2 



p+2 yi^(ex) __)r f 

-'"2ii ~ |"siS3"s2S4-'l -^2 "siS2"siS3"s2S4^'''' 2;x'( 1 — x) ( 1 — x') 



(266) 



+im\/2(a;' — x) 



Si 



S2 



XX 



{l-x){l-x') 



Whereas in the annihilation channel we calculate 

a-k\ 



-P [xtr*(-P4,P3,-(pi+P2))xs4] = X 



S3 



X' 



a' + a' 



a-k\ 1 

+ tm^Tz -cr~ 



1 - x' x'{l - x') 
x'{l-x^ 



S3SA 



Xs4 

(267) 



and 

P^bcf^'iPU -P2,Pl+P2)X 



Sli 



X 



S2 



a ■ k\ 



■a ~ a 



1 — X X 

= TlSs^-s^ - im—^ — -( 
x{l — X) 



^cr-k± . 1 i 
im^z rcr 



x{l — x) 
■^)si£l^S,,s2 , 



X. 



Sl 



(268) 



where we have introduced 



n = 



X' 1 — x' 

1 — X X 



(269) 
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We finally have 



^3132^8384^3 ■ ^4 + 5siS2'^S3S4'^SlS32^ XX'{1 — — X') 



(270) 



' x(l — x) 



''x'{l - x')' 



C Fermion and photon self energy terms 

We calculate here the fermion and photon self energy terms, arising from the second order 
commutator [rj^^^ , -^667] • 

I. We first derive the electron self energy terms. Making use of the expressions for the 
generator of the unitary transformation r/^^^ defined in eq. ( |118|) and of Hee'y, eq. (0), we 
obtain the following expression for the commutator in the electron self energy sector 



2^^PlP29p2Pl Vp2Pl9piP2/ 



'e{ptf-^^^e{pt)blb,,xt2X82 



Pi -Pi 
OjPi - pt 
Pi -pt 



^ipt)KAiXtX8, 



M2ij{Pl,P2)S 



where 



M2ij{pi,P2) = r(pi,P2,Pl -P2)'^^{P2,Pl,P2-Pl] 



(271) 



(272) 



and the momentum integration over pi,p2 is implied; 1/2 stands as the symmetry factor. The 
matrix element of the commutator between the free fermion states is 



< Pi, Sl\[n^^\ Hee^]\pi, Si > self energy 

iVpiP2 9p2Pi - Vp2Pi9piP2) ^{Pt)-^ ^ M2ii{pi,P2) 

'P2 Pi — P2 



(273) 



where the integration /p is defined in eq. (|35|). We use the expression for the generator rj through 
the coupling, namely 



1 



VpiP29p2Pl Vp2Pl9piP2 



A. 



9p 



^9p2Pl 



1P2 



PlP2 



dl 



+ 9p: 



d9 



P1P2 



•2P1 



dl 



(274) 



Change of the variables according to 



Pi 

P2 

Pi -P2 



P 

Pk 

k 



brings the integral in eq. ( |273|) to the standard form of loop integration 



/ (j]p,p—k9p^k,p Vp—k,p9p,p—k) ^(p k 
Jk 



k+ 



M2ii{p,p- k) 



(275) 



(276) 



According to eq. (|59|), the integral of the commutator [r]'^''\Hee'y] defines the difference 
between the energies (or energy corrections) dp^^ — Sp^i^. Making use of 

/■'a , 1 

P1,P2,A5'p2,P1,A 9pi,P2,>^9p2,Pl, 

a) (277) 

"''a Pi — P2 — [Pi ~ P2) 
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we have the following explicit expression: 

xT\p - k,p, -k)T\p,p - k, k) 



P 
exp 



k — {p — ky 
-2 



(278) 



A ^ 2^ 



A 



exp 



A 



p,p—k 



where the solution for the ee7-coupling constant was used. Therefore the electron energy 
correction corresponding to the first diagram, fig. (BI), is 



Sp 



lA 



(fk^dk+ e{k^ 



-e{p^ 



2(27r)3 A;^ 
xT'^p - k,p, -k)T\p,p- k, k) 



(279) 



1 



x(-i?) 



p — k — {p — k)' 

where we have introduced the regulator R, defining the cutoff condition (see main text), 



i? = exp <^ -2 



(280) 



(note that Ap^^ = Ap^p^k)- To perform the integration over k = (A;+, fc-*-) explicitly, choose the 
parametrization 



p^ 



X 



k = (xp"*", xp"*" + K"*") 



(281) 



where p = {p^,p^) is the external electron momentum. Then the terms occuring in Sp^^^ are 
rewritten in the form 



rip - k,p, -k)rip,p- k, k) = ((4^, - 4i + 2) 4 + 2m' X 



Ap^p-k =P - k - {p-ky 



\_s (x(l — x)p^ ~ xrn?) 



i^x'^ 



(282) 



Therefore the integral for the electron energy correction corresponding to the first diagram of 
fig. (0) takes the form 



p^6p 



^TT^ JO 



dx / dn 



2 



+ 1)/?^ + m'x 



2^2 



(l-x)(Ki + /(x)) 



x(-i?) 



^vr^ Jo 



dx dn. 



X 



p'^ — m' (2 



^ -2 + x 



2m' 



2 1 

+ 



4 + /(a;) \W [^-A 



(283) 



x(-i?) 



where 



/(x) = xm^ — x(l — x)p' . 



(284) 



In the last integral the principal value prescription for as x — > was introduced (see main 
text), to regularize the IR divergencies present in the longitudinal direction. 

We thus have derived the expression for the energy correction which has been used in the 
main text. 
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II. We repeat the same procedure for the photon self energy. The second order commutator 
Hee-y] gives the following expression in the photon self energy sector 

7^iVpiP2 9p2Pi ~ Vp2Pi9piP2) ' ^iPi)^i~P2)~r^ +r (^-q(^-q^x^\ (285) 
^ L yPi ~ P2 ) 



{P2 ~ Pi ) 



TrM2ij{pi,P2) 5q, -(P1-P2) 



where M2ij{pi,p2) is defined in eq. ( p72| ) and the trace acts in spin space; the integration over 
the momenta q, pi and p2 is implied. The matrix element between the free photon states reads 

< q,\\[r]'^^\Hee-y]\q,X > self energy ^ij (286) 

= — V / (^pipafepi - VP2Pi9p^p2) 0{-pt)0{pt) TrM2ij{pi,P2) 5g,-(pi-p2) , 
q^ Jpi,P2 

that can be rewritten after the change of coordinates according to 

Pi = -k 

P2 = -{k - q) 

P2-P1 = q (287) 
in the following way 

^ J^{T]k,k-ggk-g,k - r]k-q,k9k,k-q) 6'(A;+)6'(g+ - A;+) TrM2ij{k, k - q) , (288) 
where the symmetry 

V-Pl,-P2 ~ ~Vpi,P2 

9~pi,-P2 ~ 9pi,P2 (289) 
has been used. The integration of the commutator over / in the flow equation gives rise to 

iS,T. - ^,1^^-' - ^e' I '^Oik^'m,^ - k-, (290) 



X 



Tr (r(fc, k-q, q)T^{k - g, k, -q) 



exp { -2 ( ^)'| - exp |-2 ' 



This means for the photon energy correction 

5q-,,5^^ = ^^'J^^^(^ik^n^^-k^) (291) 

X Tr [r{k, k-q, q)r{k - q, k, -q)) x (-/?) , 

where the regulator R 



A, 



i? = exp|-2(^=^j I (292) 
has been introduced. Define the new set of coordinates 

(g - ky 

T = X 

q+ 

k = ((I - x)q+ , {1 ~ x)q^ + K^) 

q-k = {xq+,xq^ - K^) , (293) 



2" 
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where q = {q~^, q^) is the photon momentum. Then the terms present in (5g^;^ are 
rik, k - g, q)rik - q, k, -q) = (^rp^ {{^x - 2 + ^) + f ) 



A.-.. = q--k--{q- k)- = + $ = ^ ■ (294) 

The integral for the photon energy correction corresponding to the first diagram of fig. ([7p takes 
the form 

with 

/(x) =m2-g\(l -x) , (296) 

and the principal value prescription, denoted by '[ ]', introduced to regularize the IR divergen- 
cies. 

This is the form of the photon correction used in the main text. 
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particle number energy 

1 1 1 

► 



Figure 1: Pentadiagonal form of field theoretical Hamiltonian in 'particle number' space (the 
case of one-body Hamiltonian). Not squared region corresponds to zero matrix elements. 




Figure 2: Similarity renormalization (SR) in 'energy space' results to the band-diagonal form for 
effective Hamiltonian, where all matrix elements are sqeezed in the energy band \Ei — Ej\ < A. 



Figure 3: Similarity renormalization (SR) in 'particle number space' results to the block- 
diagonal form for effective Hamiltonian, where each block (sector) conserves the number of 
particles. 
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Figure 4: Modified similarity renormalization (MSR) of field theoretical Hamiltonians combines 
both similarity renormalization in 'particle number' and 'energy' representations. For the finite 
value of A (after the unitary transformation U{A oo, A) is performed) the matrix elements of 
'particle number changing' sectors are sqeezed in the energy band \Ei — Ej\ < A on the left hand 
side picture and are eliminated completly as A — > (that corresponds to ?7(A — > oo, A — > 0)) 
on the right hand side picture. One ends up with the block-diagonal in 'particle number space' 
renormalized effective Hamiltonian. 
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Table 1: The renormalized to the second order effective hght front QED Hamiltonian matrix in 
the Fock space. The matrix elements of the 'diagonal' (Fock state conserving) sectors exist for 
any energy differences; the 'rest' (Fock state changing) sectors are sqeezed in the energy band 
^PiPf ~ I ~ SP/ I < \ black dots correspond to the zero matrix elements to the order 
O(e^). Instantaneous diagramms are not included. 
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i = 1,2 



(7 ■ P2 — imx 



\pt] 



-a — a 



;cr - pf + imx 
\pt] 



Pi 



P3 



P2 



Pi 



Pi ki (i) 
^AAAAAA/ 



WVWW 
k2 (j) 



P2 



2 + i* i 



\pt-kt 



Pi 



P2 



Pi 



■ 



P3 



Pa 



ki (^) 



k2 U) 



P2 



1 



[Pl - pf. 



' \ -i-A \/ CA^ -i-A^ 

^PIPSA "T ^P4P2A ) [ _ J PIP3A ^ '-'P4P2A 



exp 



A2 



M2ij,X^ (xt^\{Pl,P3,Pl-P3)Xl 

X (X2'rj^(-P4, -P2, -(P1-P3))X4 



exp 



+ A^ 
A2 



M2ij,x = X^r\(pi,pi-A;i, /ci) r^(pi-A;i,p2, A;2)xi 
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Pi 



P2 



k2 U) 



P3 
— 



P4 




-el exp 



1 

X- 



A2 1 1 



A2 



\pt-pt] 



1 1 
h 



1 - exp <^ -2 



M2ij,\ = (x2'r\(pi,p2,Pl-P2)Xl 

x(X4rA(P3,P4, -(Pl-P2))X3 



A2 



e^exp 



A2 



1 

X- 



+ 



1 



2 ApjfcaA, 



1 - exp <^ -2 



A2 



M2ij,x = xt'^x{Pi,Pi-ki,ki) T{{pi-ki,p2,k2)xi 



Figure 5: The matrix elements of the renormahzed to the second order effective Hamitonian 
together with corresponding diagrams. The diagrams 2 — 5 belong to the 'diagonal' sector; the 
1,6,7 correspond to the 'rest' sector (the 6,7 diagrams are drawn schematically, namely the 
corresponding momentum change must be performed to get the real 'rest' diagrams, depicted 
in Table 1.) The photon momenta arc x+-ordered, from left to right. The similarity function 



is chosen /p,p^,A = exp( 



-A^p^^X^), where A 



—^Pf (the index 'z' denotes initial and 



'/' final states) and Ap^p^^^Pi -P2 -{P1-P2) , P ={p']_+m 
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+ 



p 



p 




+ 




Figure 6: Electron self energy: the first diagram corresponds to the commutator term [r]^'^\ H^^,^] 
in the electron self energy sector, next two diagrams arise from the normal ordering of instan- 
taneous interactions. 



k 



P 

vwww 




+ 



p 



p 




Figure 7: Photon self energy: the first diagram comes from the commutator [r]^-^\ Hee-y] in 
the photon self energy sector, the second one from the normal ordering of the instantaneous 
interaction. 
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Pi {x, k^) 



P3 {x', k'^) 



+ 



P2{l-X,-k^) P4{l-x',-k'^) 



+ 




+ 




Figure 8: The renormalized to the second order effective electron-positron interaction in the 
exchange channel; diagrams correspond to generated, instantaneous interactions and two per- 
turbative photon exchanges with respect to different time ordering. 



61 



